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I .  INTRODUCTION 


A.  THE  IONOSPHERE 

The  upper  atmosphere  is  ionized  by  solar  radiation.  As 
the  radiation  penetrates  the  atmosphere,  it  encounters  in¬ 
creasing  density  of  ionizable  material  but  is  increasingly 
attenuated  by  the  atmosphere  above,  and  these  opposing  effects 
produce  a  maximum  in  free  electron  density  at  about  300  km. 

This  ionized  region  of  the  atmosphere  is  known  as  the  iono¬ 
sphere.  The  free  electrons  in  the  ionosphere  interact  with 
radiowaves  of  frequencies  up  to  100  MHz  or  so,  and  signifi¬ 
cantly  affect  their  propagation. 

Although  the  peak  of  electron  density  occurs  at  a  height 
of  about  300  km,  subsidiary  peaks  or  "layers"  occur  at  lower 
levels  due  to  variation  with  height  of  the  atmospheric  com¬ 
position. 

The  traditional  classification  of  the  ionospheric  regions 
and  layers  as  given,  for  example,  by  DAVIES  (1965),  is: 

Height  range  (km)  Region  Layers 

50-90  D  D 

90  -130  E  E1#  E2,  Es 

130  -  up  F  Fj,  F^,  F2 

The  simplest  possible  type  of  ionized  layer  that  has  been 
derived  theoretically  is  the  so-called  Chapman  layer  which  re¬ 
sults  from  assuming: 

a)  an  atmosphere  with  only  one  type  of  gas 

b)  plane  stratification 

c)  parallel  rays  of  monochromatic,  ionizing,  solar  radiation 
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d)  an  isothermal  atmosphere „ 

It  can  then  be  shown  that  the  rate  of  production  of  ion 


pairs  when  the  angle  to  the  sun  from  the  zenith  is  x»  is  given 
at  level  c  by: 


q<X»  c)  *  q0  exp  Cl  -  C  -  sec  X  exp  (-5)]  (1) 


where  qQ  is  the  maximum  rate  of  production  of  ion  pairs  at 
level  c  *  0  when  the  sun  is  overhead,  and 


(2) 


where  h  is  the  height  above  the  ground,  hQ  is  the  height  of 
maximum  ion  production  and  H  is  the  scale,  height  of  the  iono¬ 
sphere  given  by 


H  - 


kT0 

Mg 


(3) 


where  k  is  Boltzmann's  constant,  Tq  is  the  uniform  temperature 
of  the  atmosphere,  M  is  the  mean  molecular  mass  of  the  atmospheric 
gas,  and  g  is  the  constant  acceleration  of  gravity. 

Assuming  that  the  removal  of  electrons  is  due  to  a  recom¬ 
bination  process  described  by 

jjr  -  q  -  oNJ  (4) 

where  N  is  the  electron  number  density,  t  is  time,  and  a  is  the 
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recombination  coefficient;  and  also  assuming  that  equilibrium 
has  been  reached,  there  results: 

N(C)  =  NQ  exp  j  [1  -  ?  -  sec  x  exp  (-?)]  (5) 

where  N  is  the  maximum  electron  number  density  at  level  ?  =  0 
o 

when  the  sun  is  overhead.  This  expression  is  known  as  the 
Chapman  law.  In  its  derivation,  the  recombination  coefficient 
is  taken  to  be  independent  of  height.  In  view  of  this  simpli¬ 
fying  assumption,  and  the  four  previous  ones,  it  is  surprising 
that  the  B  and  F^  layers  do  in  fact  approximately  follow  the 
Chapman  law.  In  Figure  1,  the  expression  for  N (5)  is  plotted 
against  C  for  several  values  of  the  solar  zenith  angle  x* 
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B.  r,U)IO  WAVES  IN  THE  IONOSPHERE 

The  existence  of  the  ionosphere  was  postulated  by 
Kennelly  and  Heaviside  in  1902  to  account  for  the  success 
of  Marconi's  trans-Atlantic  radio  experiment.  However  its 
existence  was  not  confirmeJ  beyond  reasonable  doubt  until 
1925  when  in  a  single  year  six  independent  methods  of  measur¬ 
ing  the  height  of  the  ionosphere  were  published  and  (or)  used. 
Most  of  the  knowledge  now  existing  about  the  ionosphere  has 
been  and  is  being  obtained  by  means  of  radio  techniques.  It 
is  appropriate,  therefore,  to  review  the  general  aspects  of 
radio  waves  propagating  in  the  ionosphere.  Perhaps  the  most 
complete  treatment  of  this  topic  is  that  by  3UDDEN  (1961). 
Generalized  treatments  of  radio  waves  in  magneto-plasmas  are 
also  available  such  as  that  by  STIX  (1962). 

In  the  presence  of  the  earth's  magnetic  field,  the  iono¬ 
spheric  magneto-plasma  acts  as  a  bi-refringent  medium,  decom¬ 
posing  a  plane  polarized  wave  into  "ordinary"  and  "extra¬ 
ordinary"  waves  and  changing  the  direction  of  energy  flow 
away  from  the  direction  of  phase  propagation. 

For  the  present  purpose  it  will  suffice  to  limit  con¬ 
sideration  to  the  propagation  of  radio  waves  of  the  order 
of  a  few  megahertz,  and  it  will  not  be  necessary  to  take 
into  account  the  motion  of  the  neutral  particles,  the  pres¬ 
sure  gradients  of  the  atmospheric  constituents,  or  the  fric¬ 
tional  effects  of  collisions.  With  these  simplifications  and 
the  assumption  of  small  plane  wave  disturbances  (entirely 


-  f ,-  _1*@\  ¥  ;  ~r-JKr  < 
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equivalent  to  cold-plasma,  collisionless  wave- theory)  it  can 
be  shown  that  the  dispersion  equation  of  the  characteristic 
waves,  i.e.,  the  waves  which  can  propagate  in  the  medium, 
is  given  by: 


n2 


_ 2X(i-x) _ _ 

2  (l-X)-y2sin20±/5f  4sin'*0+4  (1-X)  2y2cos*6 


(6) 


which  is  known  as  the  Appleton-Hartree  equation.  Here  n  is 
the  refractive  index  of  the  ionosphere  for  the  radio  wave; 

8  is  the  angle  between  the  earth's  magnetic  field  and  the 
direction  of  the  wave-normal; 


(7) 


where  w  is  the  angular  frequency  of  the  radio  wave,  is 
the  local  angular  plasma  frequency  given  by 


2  Ne2 
^p  =  eom 


(8) 


where  eQ  is  the  permittivity  of  free  space,  e  is  the  charge 
of  an  electron(  m  its  mass;  and 


y  » 


(9) 


where  o>H 


where  13 

o 


is  the  electron  angular  gyro frequency  given  by 
eB 

u)„  =  — -  (10) 

h  m 

is  the  earth's  magnetic  field,  taken  to  be  constant. 
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In  general,  the  Apple ton-Har tree  dispersion  relation 
yields  four  possible  values  for  the  refractive  index  n,  once 
a  particular  set  of  values  for  X,  Y,  and  0  has  been  chosen. 

These  occur  in  pairs  which  have  the  same  magnitude  but  opposite 
sign.  Each  of  these  two  pairs  describes  different  types  of 
waves,  whereas  each  component  of  a  given  pair  represents  the 
same  type  of  wave  but  with  opposite  directions  of  propagation. 

In  addition  to  the  phase-velocity  described  by  the  Appleton- 
Hartree  dispersion  relation,  it  is  useful  to  consider  the  polar¬ 
ization  of  the  characteristic  waves.  The  polarization  of  a 
plane  wave  travelling  in  the  z-direction  in  the  ionosphere  is 
defined  by  the  ratio  R  of  the  transverse  electric  field  com¬ 
ponents 


R  = 


(ID 


when  the  coordinate  system  is  chosen  so  that  Bq  lies  in  the 
yz-plane.  Using  the  parameters  of  the  dispersion  relation, 
the  polarization  becomes: 


R  -  - 


2Y  cos  6 


sxn' 


0  —  iZ  -J  v  .  , 

+^[T  ,-ySi  +  4  cos 


sin4  0 


(12) 


The  reflection  condition  for  the  case  of  a  radio  wave 
Vertically  incident  upon  the  ionosphere  is  obtained  by  setting 
n  *=  0  in  the  dispersion  relation  and  solving  for  X, 
with  4  sign:  X  =  1 
with  -  sign:  X  =*  1  -  Y 


or 


X  *  1  +  Y 


(13) 

(14a) 

(14b) 
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Equation  (13)  shows  that  for  one  of  the  characteristic  waves, 
the  wave  is  reflected  as  in  the  absence  of  a  background  mag¬ 
netic  field,  i.e.,  Y  «=  0.  This  wave  is  called  the  ordinary 
wave.  Equations  (14)  show  that  the  height  of  reflection  of 
the  other  wave  (colled  extraordinary)  is  dependent  on  the 
magnitude,  but  not  on  the  direction  of  the  earth's  magnetic 
field.  For  waves  of  frequency  greater  than  the  gyrofrequency 
the  height  of  reflection  is  given  by  X  =  1  -  Y  whereas  for 
frequencies  below  the  gyrofrequency  the  condition  is  X  =  1  +  Y. 


C.  THE  IONOSONDE 


The  most  widely  used  instrument  for  investigating  the 
ionosphere  by  radio  is  the  ionosonde.  This  is  essentially 
a  pulsed  radar  device  in  which  the  .probing  frequency  is 
swept  from  about  2  MHz  up  to  20  MHz. 

The  equipment  is  designed  to  measure  directly  the  time 
taken  for  a  pulse  of  radio  waves  to  travel  up  to  the  iono¬ 
sphere  from  the  ground  (or  down  to  it  from  a  satellite)  and 
back,  as  a  function  of  frequency.  This  time  is  used  to  define 
the  radar  range  to  the  reflection  ight.  Assuming  group- 
propagation  at  the  speed  of  light  it  defines  h',  the  so-called 
"virtual  height"  of  the  ionosphere,  given  by 

h'  =  0.5  ct  (15) 

where  c  is  the  speed  of  light  in  free  space  and  t  is  the  echo 
delay  time.  The  resulting  photographic  record  of  virtual 
height  versus  frequency  is  known  as  an  ionogram. 

Although  the  radar  pulse  travels  progressively  slower 
as  it  penetrates  the  ionosphere,  it  is  generally  possible  to 
compute  from  an  ionogram,  the  electron  number  density  N(h)  as 
a  function  of "true-height"  h  and  so  provide  an  electron  density 
height-profile  of  the  ionosphere. 

The  relation  between  true  height  h  and  virtual  height  h' 
is  as  follows 

h 

h ’  *  /  ji'dh  (16) 

o 


s 


where  the  group  refractive  index  y  *  is  given  by  i 

i 

y'  =  g~  (nto)  (17)  j 

j 

The  inversion  of  equation  (16)  is  complicated  by  the  \ 

\ 

fact  that  \i' i  defined  through  equations  (17) ,  (6 ) ,  (7)  #  and  ^ 

£ 

(8),  is  an  unknown  function  of  h..  The  solution  may  be  effected  ] 

in  either  of  two  ways,  as  classified  by  THOMAS  (1959):  j 

{a)  Model  or  comparison  methods  in  which  the  layer  is 
assumed  to  have  an  N  (h)  profile  such  that  substi¬ 
tution  in  equation  (16)  gives  an  integral  which  can 
be  evaluated  analytically.  The  calculated  h'  (f) 
curve  is  then  compared  with  the  observed  one  and  the 
parameters  of  the  assumed  distribution  obtained. 

(b)  Integral  equation  methods  which  may  be  divided  into 
two  categories:  in  the  first,  the  equation  is  in¬ 
verted  directly  to  give  an  analytic  solution;  and 
in  the  second,  the  integral  is  replaced  by  a  dis¬ 
crete  sum  over  a  number  of  sufficiently  small  frequency- 
intervals  . 

Model  or  comparison  methods  have  the  basic  disadvantage 
that  the  shape  of  the  layer  is  assumed  a  priori  and,  therefore, 
they  need  not  be  discussed  any  further. 

Both  integral  equation  methods  suffer  from  the  inherent  instru¬ 
mental  limitation  that  there  are  no  observations  of  h' (f)  below  some 


. .  , 
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value  of  the  probing  frequency.  Direct  inversion  methods  are 
unwieldly  except  when  the  magnetic  field  of  the  earth  is 
neglected  and  even  when  it  is  not,  they  apply  exactly  only  on 
the  magnetic  equator  (6  *=  90°).  A  second  fundamental  limita¬ 
tion  of  all  integral  equation  methods  is  that  it  is  necessary 
to  assume  that  the  plasma  density  is  a  monofconically  increas¬ 
ing  function  of  height.  A  more  recent,  simpler  and  faster 
integral  equation  method  due  to  TITHERIDGE  (1967)  has  been 
adopted  here  and  is  described  in  the  next  section. 


D.  THE  OVERLAPPING  POLYNOMIAL  METHOD  OF  ANALYSIS  OF  IONOGRAMS 
In  this  method,  the  true  heights  of  reflection  at  a 
series  of  frequencies  f2»  f ^ ,  ...»  f  are  calculated  from 
the  measured  virtual  heights  at  these  same  frequencies.  When 
computing  the  real  height  h .  at  the  frequency  >  the  real 
height  curve  between  the  frequencies  f^-2  and  rs  repre¬ 

sented  by  a  five- term  polynomial  in  the  plasma  frequency  fN- 
The  polynomial  is  required  to  pass  through  the  real  heights 
^i-2  *  ^i-1  *  f±  *  and  to  give  the  observed  virtual  height 
at  the  frequencies  and  These  conditions  define 

a  set  of  five  polynomial  coefficients  at  each  frequency ,  giving 
the  value  of  h^  in  terms  of  hi-2 ,  hi-i»  h'i-i»  h'i»  and  hW 
The  polynomial  used  in  calculating  h^  extends  only  down 
to  the  height  h^_2>  The  contribution  of  the  ionization  below 
this  height  to  the  values  of  h'^_^,  h'^,  h'i+^  must  therefore 
be  allowed  for  separately.  This  is  done  by  the  method  used  in 
the  normal  lamination  analysis  of  ionograms ,  i.c.  ,  by  the  re¬ 
placement  of  integral  (16)  by  a  discrete  sum  over  a  number  of 
sufficiently  small  frequency  intervals. 
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E.  TRAVELLING  IONOSPHERIC  DISTURBANCES 

Iso- ionic  contours  are  height-contours  of  constant  elec¬ 
tron  density  in  the  ionosphere.  Perturbations  of  the  normal 
equilibrium  conditions  create  ripples  of  horizontal  gradients 
in  these  contours.  Anomalous  reflections  then  appear  on  iono- 
grams  and  records  of  other  radio  wave  observing  stations  de¬ 
pending  upon  ionospheric  reflection.  Spaced  receiving  sites 
often  exhibit  time  delays  between  them  in  the  occurrence  of 
similar  phases  of  perturbations ,  thereby  demonstrating  hori¬ 
zontal  propagation  of  them. 

The  most  successful  interpretation  of  these  observations 
is  that  which  attributes  travelling  ionospheric  disturbances 
(T.I.D.'s)  to  corresponding  perturbations  of  the  neutral  gas 
associated  with  the  passage  of  internal  atmospheric  gravity 
waves.  This  interpretation  is  consistent  with  T.I.D.'s  in 
so  many  respects  that  there  is  little  doubt  about  its  basic 
correctness#  HOOKE  (1968). 
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F.  INTERNAL  ATMOSPHERIC  GRAVITY  WAVES 

The  theory  of  internal  gravity  waves  has  been  studied 
extensively  by  HINES  (1960).  The  analysis  is  based  on  an 
idealized  model  in  which  the  atmosphere  is  taken  to  be 
stationary  in  the  absence  of  such  waves,  and  to  be  uniform 
in  both  temperature  and  composition.  Superimposed  wave 
motions  are  assumed  to  have  only  perturbation  magnitude,  in 
that  terms  of  second  and  higher  orders  are  neglected  in  the 
equations  governing  the  motion.  Changes  are  assumed  to  occur 
adiabatically  and  the  gravitational  field  is  taken  to  be  con¬ 
stant  both  in  direction  and  magnitude. 

Under  these  assumptions ,  it  can  be  shown  that  atmo¬ 
spheric  oscillations  are  governed  by  the  following  dispersion 
relation : 

w*  -  u>2C2  (k£  +  k2)  +  iw2y.gkz  +  (y-l)g2k^  =  0  (18) 

where  w  is  the  angular  frequency  of  the  wave,  y  is  the  specific 
heats  ratio  of  air,  g  is  the  acceleration  of  gravity,  k^  and 
kz  are  the  horizontal  and  vertical  components  of  the  phase 
propagation  vector  Jc  (or  angular  wc  *e  number) ,  and  C  is  the 

S 

speed  of  sound  given  by 

C2  =  yHg  (19) 

where  K  is,  again,  the  scale  height  of  the  atmosphere  as 
defined  in  equation  (3) .  Hines  confined  his  attention 
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to  the  case  that  k^  is  purely  real  so  that  equation  (18) 

reveals  two  possibilities:  either  k  is  purely  imaginary  or 

z 


+  V2H 


where  kzr  is  purely  real. 

The  first  alternative  corresponds  to  horizontally  propa¬ 
gating  "surface"  waves  having  vertical  attenuation.  The 
second  alternative  corresponds  to  "internal"  waves  which  pos¬ 
sess  a  vertical  component  of  propagation.  Substituting  equa¬ 
tion  (20)  into  (18)  results  in  the  dispersion  relation  for 
internal  waves  which  may  be  cast  in  the  following  form: 


w* 

1  ~  -f  , 

_ M±  + _ *„  „„ 

1  -  5^  1  -  5T 


n2  =  1 
zr 


where 


»  JT_-  Dj 

g 


•4-8 


"h  -  — 


cskh 


ck 

s  zr 


zr  w 


Since  Y  is  less  than  2  it  follows  that  u>  <  w  -  From 

g  a 

equation  (21)  it  can  be  shown  that  there  exists  a  gap  in  the 
frequency  spectrum 


tog  <  to  <  W. 


(26) 


15 


in  which  no  internal  waves  can  be  propagated.  This  indi¬ 
cates  that  two  distinct  sequences  of  internal  waves  can 
occur.  The  one  at  high  frequencies,  w  >  tn  (T  <  T  ) ,  and 
the  other  at  low  frequencies,  w  <  <ug  (T  >  T^)  where  T  =  2ir/w 
and  Tg  =  2Tr/Wg.  The  two  sequences  have  been  termed  by  Hines 
"internal  acoustic"  waves  and  "internal  gravity"  waves. 

In  Figure  2,  contours  of  constant  period  in  the  n^  -  nzr 
domain  are  shown  for  the  internal  gravity  wave  mode.  From 
these,  the  directions  of  phase  propagation  and  energy  flow 
can  be  readily  determined,  STIX  (1962) .  The  direction  of 
phase  propagation  is  that  of  the  radius  vector  (n^,  nzr) 
drawn  from  the  origin  to  any  particular  point  in  the  n^  vs. 
n  domain.  The  direction  of  energy  propagation  is  that  of 
the  normal  to  the  corresponding  constant  period  contour, 
drawn  at  the  tip  of  the  radius  vector  (n^,  nzr) .  It  is  seen 
that  energy  m&y  propagate  in  directions  radically  different 
from  the  corresponding  phase  normals,  by  as  much  as  90°  when 
the  asymptotes  of  the  period-contours  are  approached,  and  for 
the  longer  periods  the  flow  tends  to  be  nearly  horizontal. 
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G.  TRAVELLING  IONOSPHERIC  DISTURBANCES  AND  GRAVITY  WAVES 

In  a  perturbation  treatment  of  the  ionospheric  irregu¬ 
larities  caused  by  internal  atmospheric  gravity  waves,  HOOKE 
(1968!  arrived  at  several  important  conclusions.  In  the  F2- 
region,  the  main  effect  of  the  gravity  wave  is  to  impart  the 
component  of  motion  of  the  neutral  gas  which  is  parallel  to 
the  magnetic  field,  to  the  ionized  component.  Up  to  the  Fl- 
ledge,  the  effect  of  gravity  waves  on  the  rate  of  chemical 
processes  is  important  and  at  any  point,  gravity  waves  affect 
the  rate  of  photoionization  by  altering  both  the  neutral  gas 
number  density  and  the  ionizing  radiation  flux  at  that  point. 

At  heights  at  or  above  the  height  of  the  F2-peak,  the 
amplitude  of  the  perturbation  in  electron  number  density  N’  at 
(x,  y,  z,  t)  is  given  by 

|N'|-  No(z)Ub(zQ)sin  I  exp[kzi  (z-zo)]w_l 

(27) 

assuming  that  it  is  caused  by  an  internal  atmospheric  gravity 
wave  producing  neutral  gas  velocities  along  the  earth’s  magnetic 
field  given  by 

«b  *  Ub^)eKpCkzi{z~2o)JexP  “  kxx  ”  kyY  -  \rz) 

(28) 


where 


V  s  _±_ 

Kzi  2H 


(29) 


II .  THEORETICAL  WORK 

A,  GRAVITY  WAVE  RESONANCE  CONDITION  jc  •  $  =  0 

o 

There  is  ample  evidence  that  travelling  ionospheric 
disturbances  (T.I.D.'s)  are  caused  by  the  passage  of  inter¬ 
nal  atmospheric  gravity  waves,  HINES  (I960). 

The  nature  of  the  interaction  of  gravity  waves  with 
the  ionosphere  has  been  studied  in  detail  by  HOOKE  (1968) . 
Previously,  HINES  (1955)  had  sought  to  explain  T>I;D.’s 
as  caused  by  the  resonant  response  of  the  ionosphere  and 
found  two  conditions:  3c  *  k  =  0  and  k  •  Bq  =  0,  where  3c 
is  the  gravity-wave  complex  propagation  vector  and  3q  is 
the  magnetic  field  of  the  earth.  Since  k  is  in  general 
complex,  ic  •  Jc  =  0  does  not  necessarily  imply  that  ic  =  ($. 

The  condition  ic  *  BQ  «  0  means  that  3c  must  lie  in  a  plane 
perpendicular  to  Bq. 

Neither  author  analyzed  the  case  k  *  Bq  =  0  which  is 
the  subject  of  the  present  section.  Hines's  gravity  wave 
dispersion  relation  is  generalized  so  as  to  include  hori¬ 
zontal  attenuation  or  growth.  It  is  then  found  that  the 
frequency  gap  between  the  internal  "acoustic”  and  "gravity" 
modes  persists.  However,  the  internal  gravity  wave  mode  no 
longer  encompasses  all  frequencies  under  the  Vaisala-Brunt 

frequency,  but  becomes  a  band  limited  by  w  sin  A  and  »u  , 

y  y 

where  A  is  the  magnitude  of  the  magnetic  dip  angle. 
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1.  Mathematical  Development 

The  basic  assumptions  are  those  of  HINES  (1955).  The 
atmosphere  is  taken  to  be  an  isothermal,  stationary  medium, 
permeated  by  a  uniform  magnetic  field  $o  and  subject  to  a 
constant  gravitational  acceleration  g. 

The  analysis  is  based  on  three  fluid  equations :  the 
conservation  of  mass,  the  adiabatic  equation,  and  the  equation 
of  motion;  and  three  electromagnetic  equations:  the  two 
Maxwell's  equations  and  a  tensorial  Ohm's  law. 

In  the  equation  of  conservation  of  mass,  the  source 
term  is  set  equal  to  zero;  in  the  equation  of  motion  the 
effect  of  gravity  is  included  but  the  electric  and  magnetic 
forces  are  set  equal  to  zero. 

In  Ampere's  law,  the  displacement  current  is  neglected, 
and  in  the  conductivity  tensor  a,  the  fact  that  the  longi¬ 
tudinal  conductivity  oQ  is  several  orders  of  magnitude  larger 
than  the  transverse  conductivities  will  be  used. 

HINES  (1955)  showed  that  the  foregoing  conditions  amount 
to  assuming  that  hydrodynamic  oscillations  of  the  ionosphere 
dominate  the  basic  disturbance,  though  secondary  electromag¬ 
netic  effects  may  play  an  important  part  in  the  observations. 
The  medium  considered  is  the  simplest  one  but  still  retains 
the  essential  features  of  the  ionosphere. 

Under  the  conditions  listed  above,  the  pertinent  equations 

are : 

ff  +  7  •  (pv)  =  0  (30) 

&  -  0 


{3 1) 
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p  +  pg  (32) 

V  X  2  =  -  ||  (33) 

7  x  S  =  yo?  (34) 

3  =  a  •  (E  +  v  x  Bq)  (35/ 


where  p  is  pressure,  p  is  density,  v  is  the  velocity  of  a 
fluid  element,  y  is  the  specific  heats  ratio  of  the  fluid, 

2  and  il  are  the  electric  and  magnetic  fields,  is  the 
magnetic  permeability  of  free  space,  and  J  is  the  current 
density. 

The  geometry  adopted  here  is  shown  in  Figure  3.  The  z~ 
axis  is  taken  in  the  vertical  direction  with  g  =  -  gz, 
is  a  constant  magnetic  field  in  the  yz-plane  and  jc  is  a  propa¬ 
gation  vector  in  an  arbitrary  direction. 

Assuming  now  that  small  changes  of  the  form  exp  [i(tot-k*r)3 
take  place  in  all  variables  but  that  the  pressure  and  density 
also  vary  as  exp  (-  z/H) ,  where  H  is  the  scale  height  of  the 
atmosphere,  the  fluid  equations  yield 

[w 2 1  -  C2&  -  iicg  -  i  (y-l)gic]  •  v  =  0  (36) 

5 

where  f  is  the  unit  tensor,  and,  again 

C2  -  ygH  (19) 

defines  the  speed  of  sound  in  the  medium.  Equation  (36)  can 


w*r 
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also  be  writtu  :i  as 

w*  -a)2  C2(k2+k2+k2)  +  (Y-l)g2(k2+k2)  +  iYgu>2kz  =  0  (37) 

which  is,  of  course,  Hines's  dispersion  relation.  Applying 
the  same  analysis  to  the  electromagnetic  equations  it  follows 
that 


(k2I  -  iclc  +  iojyo0)  «  S  =  iwvio0oI  *  Sq  •  v  (38) 

HINES  (1955)  showed  that  hydromagnetic  resonance  occurs  for 
either  of  the  case J 


k  •  k  =  0 


(39a) 


or 


it  -  S  =o 

o 


(39b) 


which  are  obtained  by  setting  the  determinant  of  the  left 
hand  side  of  equation  (38)  equal  to  zero  and  neglecting  the 
transverse  conductivities . 

Case  a)  was  discussed  by  HINES  (1955)  and  need  not  be 
repeated  here. 

In  order  to  simplify  the  algebra,  the  following  definitions 
airs  recalled 


u 


w 


2  =  iizils: 

3  c! 

2  *  12 
a  4H 


(22) 

(23) 


■i><TiwnMar  ffr>T*%i-fr 
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also 


n  = 


C_s^ 

to 


(40) 


so  that  equation  (37)  becomes 


1  - 


+ 


i 


0 


(41) 


Calculations  become  even  simpler  using  a  new  coordinate  system 
(£,  Tl/  C)  obtained  by  a  counterclockwise  rotation  of  the  (x,y,z) 
system  about  the  x-axis  through  the  dip  angle,  A.  The  £-  and 
x-axes  are  identical.  The  n~axis  lies  along  the  magnetic  field 
and  the  ?-axis  completes  the  right  handed  rectangular  sys¬ 
tem.  Equations  (39b)  and  (41)  then  become,  respectively, 


nn  =  0 


(42) 


and 


or  01 

1  -  (n£+n2)  +  -f-  (n£+n2sin2A)  +  i  2  —  n  cos  A  =  0  (43) 


Assuming  that  n^  is  real,  then  n^  must  be  complex,  thus 


”c  =  v  +  1  V 


(44) 


where 


and 


“a 

—  cos  A 

_  _  « _ 

£1  0)? 


(45) 


1  -  —2-  sin2  A 


n2  + 


0)2 

< 

1 - y 

O)2 


or 


1  -  — §•  sin2  A  -  —  cos2  A 


n2 


0)' 


0)  1 


1  -  -§•  sin2  A 
0)2 


w*  2 

(1  -  -§■  sin2  A) 


0)* 


(46) 
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Note,  that  equations  (45)  and  (46)  reduce  exactly  to  those  of  inter¬ 
nal  gravity  waves  when  A  =  0°  except  that  equation  (46)  would  apply 
only  to  the  vertical  plane  perpendicular  to  the  magnetic  field. 
Careful  inspection  of  equation  (46)  reveals  again,  in  close 
analogy  with  Hines's  internal  waves,  that  two  modes  of  waves 
may  be  recognized:  "gravity"  and  "acoustic".  The  following 
tabulation  delineates  the  situation: 

Comments 
no  solution 
gravity  mode 
A  no  solution 

acoustic  mode 

Figures  4  and  5  show  the  variation  of  as  a  function 
of  magnetic  angle  and  period  for  selected  periods  and  angles, 
respectively. 

Figures  6  through  8  show  plots  of  n  vs.  nF  for  selected 
periods  and  for  angles  of  15° ,  45° ,  and  75° . 


mwu«n.M  n-u^ u  i.i  •wmmmmmmun* 
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2.  Results 

This  simplified  analysis  shows  the  feasibility  of  the 

£  .  S  =0.  The  main  feature  of 
resonant  gravity  wave  mode  k  bq 

this  mode  is  its  confinement  to  a  band  »g  sin  A  <  «  <  “g 
which  rules  out  periods  greater  than  .g  Csc  A  from  magnetic 
latitudes  greater  than  A.  Away  from  the  magnetic  equator, 
the  propagation  band  is  rather  narrow,  e.g.,  at  A  -  30-  rt 

is  only  T <  T  <  2Tg. 

The  indices  of  refraction  n^  and  n^r  are  entirely 
gous  to  those  given  by  HIKES  (1960)  so  that  the  phase  and 
group  velocities  still  bear  the  same  relationship:  for  the 
longer  periods,  the  energy  flow  is  nearly  horizontal,  radi¬ 
cally  different  from  the  corresponding  phase  flow.^  ^ 

As  opposed  to  resonance  mode  t  •  £  -  0,  mode  k  *  Bo  ~  0 
exhibits  both  vertical  and  horizontal  phase  shifts  and  atten- 

uation  or  growth. 


i 


Mean %, ahjhdMiS«l a  -  ■  •  ■- 


B.  IONOSPHERIC  PREDICTIVE  FUNCTION 


HOOKE  (1968)  has  shown  that  in  the  F2-region  the  ampli¬ 
tude  of  the  electron  number  density  variation  produced  by  an 
internal  atmospheric  wave  is  described  by  equations  (27)  through 
(29),  repeated  here: 

|N' |  =  No(z)Ub(zo)sin  I  exp  [kzi  (z  -  ~  +  k 

(27) 

ub  =  ub  (2o)  exp  [kzi(z  "  2on  exp  1  {u)t  "  V  "  kyy  “  kzrz) 

(28) 


The  full  electron  density  variation  is  given  by 

< 


N' 


In’*  j  exp  i 


k,  /sin  I 

ut  -  V  -  V  -  kzr  z  +  2  -  tan~‘  r~m — 


i_  +  k 

3z  zi 


(47) 


In  order  to  obtain  an  expression  for  isoionic  contours, 
the  height  variation  of  contours  of  constant  N,  it  is  neces¬ 
sary  to  set 

Nq(z)  +  Re  (N*)  =  constant  (48) 
and  solve  for  z.  This  proves  to  be  a  cumbersome  inversion 
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which  is  carried  through  only  because  one  may  wish  to  compare 
the  theoretical  expressions  with  experimental  data  which  have 
commonly  been  presented  in  the  form  of  iso-ionic  contours. 

There  is  no  fundamental  reason  to  choose  iso-ionic  contours 
over  iso-height  contours,  the  variation  of  N  at  constant  height, 
nor  greater  difficulty  presenting  the  data  under  either  scheme. 
Hence,  as  far  as  T.I.D.  work  is  concerned,  it  is  much  more 
satisfactory  to  present  experimental  data  as  iso-height  con¬ 
tours  so  that  they  can  be  readily  compared  to  the  available 
theory. 

The  derivation  of  the  gravity  wave  mode  assumed  constancy 
in  the  value  of  the  scale  height  H.  The  fact  is  that  in  the 
region  of  interest,  say  200  to  300  km,  H  varies  from  48  to 
60  km  so  that  the  theory  is  imperfect.  If,  however,  the  iono¬ 
sphere  is  approximated  by  successive  horizontal  layers  10  km 
thick# then  H  remains  constant  within  2.5%  in  each  slab,  and 
the  theory  may  be  reasonably  applied  to  each  slab. 

Equation  (27)  suggests  the  possibility  of  predicting  the 
electron  density  variation  at  one  level  when  the  variation  is 
known  at  some  other  level.  For  levels  z-j^  and  Z2 /  the  ratio  of 
the  variation  at  frequency  w  is  given  by 


where  k^r  encompasses  the  frequency  dependence. 


(49) 
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As  it  stands,  equation  (49)  may  be  thought  of  as  a  pre¬ 
dictive  function  for  an  ionospheric  layer  of  thickness 
z2  -  and  centered  about  (z^  +  z2)/2.  For  an  input,  |N'(z^)|, 
|N' (z2) |  is  predicted.  Therefore  we  may  input  an  observed 
Jn* (z1) | ,  predict  J N *  (z2) |  as  well  as  observe  Jn'(z2)|,  and 
compare  the  predicted  and  observed  values  of  | N *  (z2) |  .  This 
can  be  done  in  a  revealing  manner  by  Fourier  decomposing  the 
predicted  and  observed  functions  and  examining  their  ratio  at 
each  harmonic  frequency.  To  the  extent  that  confidence  in  the 
theory  can  be  established  in  this  way,  predictions  can  be  made 
to  heights  above  (or  below)  available  observations. 

With  this  perspective,  and  by  analogy,  Prof.  Morgan  has 
suggested  the  possibility  of  developing  a  horizontal  predictive 
function  in  which  the  observation  points  are  at  the  same  height 
but  are  horizontally  displaced. 


Ill .  EXPERIMENTAL  WORK 


A.  DARTMOUTH  IONOSONDE  NETWORK 

The  Dartmouth  ionosonde  network  consists  basically  of 
three  vertical-incidence  ionosondes  operating  independently 
at  locations  roughly  determining  an  equilateral  triangle  150  km 
on  a  side,  as  shown  in  Figure  9. 

Typical  operating  characteristics  of  these  sounders  are 
as  follows:  one  frequency-sweep  every  2  minutes,  100  pulses 
per  second,  Gaussian  pulse  envelope,  and  logarithmic  frequency 
sweep  from  about  1.6  to  20  MHz  in  approximately  20  seconds. 

The  ionosondes  were  designed  and  built  at  the  Dartmouth 
Radiophysics  Laboratory  and  a  description  of  them  will  be 
published  by  Prof.  Morgan  and  Mr.  Pratt  >o  they  will  not  be 
dealt  with  further  here.  An  outstanding  characteristic  of 
them  is  their  ability  to  run  one-week  between  film  changes 
and  their  ability  to  run  months  at  a  time  without  any  other 
attention. 

The  ionograms  are  produced  in  motion-picture  format  and 
each  weekly  film  is  inspected  by  viewing  it  as  an  accelerated 
motion  picture.  If  a  possible  T.I.D.  is  detected,  the  cor¬ 
responding  ionograms  are  digitized  by  projection  and  tracing 
on  a  Grafacon  1010A  tablet  in  conjunction  with  a  Spiras-65 
computer.  The  digital  data  are  transferred  into  paper  tape 
which  is  then  processed  by  the  Dar  tmouth  Time  Sharing  System 
to  convert  virtual  height  to  true  height  by  the  overlapping 


28 


polynomial  method  previously  described.  The  data  are  then 
plotted  by  a  Complot  DP-1-1  digital  plotter,  again  in  con¬ 
junction  with  the  Spiras-65  computer.  If  the  resulting  plot 
confirms  the  existence  of  a  T.I.D.  of  interest,  the  data  are 
further  processed  for  finer  detail  and  computation  of  speed, 
direction,  power  spectra,  and  comparison  with  the  ionospheric 
predictive  function  associated  with  the  T.I.D. 
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1.  Hardware 

In  what  follows#  a  very  brief  description  of  the  data 
processing  equipment  and  interfacing  will  be  given. 

The  ionogram  film  projector  is  a  Vanguard  M-16C  "motion 
analyzer"  which  is  a  movie  projector  capable  of  a  wide  range 
of  speed,  both  forward  and  backward  with  the  added  facility 
of  frame-by-frame  projection.  Two  front-surfaced  mirrors  are 
used  to  project  the  image  onto  the  translucent  screen  of  the 
Graf aeon  1010A  digitizer  from  beneath.  By  means  of  a  spring- 
loaded  tapered  pin  in  the  sprocket  holes,  the  projector  pro¬ 
duces  a  high  degree  of  registration  from  frame  to  frame. 

The  Grafacon  digitizer  consists  essentially  of  a  tablet 
with  a  mesh  of  imbedded  wires  and  a  stylus  capable  of  recog¬ 
nizing  its  position  on  the  tablet.  The  xy-coordinates  of  the 
stylus  position  are  given  by  two  binary  numbers  of  10  bits 
each,  picked  up  from  the  imbedded  wire  mesh. 

The  interfacing  of  the  digitizer  and  the  Spiras-65  computer 
has  involved  the  building  and  interconnection  of  twenty-three 
level- changing  circuits  for  the  twenty  data  bits  and  three 
control  bits;  additional  interconnection  within  the  computer; 
and  adjustment  of  the  power  supply  circuit  of  the  extension 
unit  of  the  computer. 

The  Complot  DP-1-1  plotter  is  a  conventional  digital 
plotter.  The  interfacing  of  the  plotter  and  the  Spiras-65 
computer  has  required  the  modification  of  the  internal  circuits 
of  the  former  and  the  addition  of  circuits  to  the  latter. 
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The  input  circuits  of  the  plotter  were  altered  to  accept 
the  computer  pulses.  Four  cards  were  added  to  the  computer: 
two  input/output  gating#  one  input/output  controller#  and  one 
special-purpose  card.  The  first  three  were  standard  cards 
while  the  special-purpose  card  was  designed  and  built  at 
Dartmouth.  This  card  allows  enough  time  fcr  the  plotter  to 
perform  its  functions. 

The  Spiras-65  computer  is  a  small  multipurpose  computer 
with  an  8k  memory  attended  by  a  model  33  ASR  teletype  unit. 

All  of  the  interfacing  circuits#  modifications,  and 
wiring  lists#  are  appended  to  the  "Small  Computer  Interfacing" 
volume  at  the  Radiophysics  Laboratory. 
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2.  Software 

The  main  data  digitizing  program  is  named  "Ionogram 
Scaling".  It  is  written  in  Spiras-65  machine  language  be¬ 
cause  only  4k  of  memory  was  available  until  recently.  This 
program  simply  allows  the  operator  to  scale  ionograms  in  an 
automatic  fashion.  Its  main  features  are  visual  time-check 
facility,  calibration  for  both  scale  and  orthogonality  in  the 
ionograms,  determination  of  the  critical  frequency,  automatic 
recording  of  virtual  heights  at  a  preset  number  of  frequencies, 
and  memory  dump  of  coded  data  on  paper  tape.  Two  slight  vari¬ 
ations  of  this  program  have  also  been  written:  "Critical 
Frequencies"  and  "2-Virtual  Heights/Ionogram" .  The  former 
records  on  paper  tape  only  the  critical  frequency  of  each 
ionogram  while  the  latter  records  any  two  virtual  heights  per 
ionogram. 

The  plotting  program  is  called  "Plot"  and  is  also  written 
in  Spiras-65  machine  language.  It  accepts  coded  data  in  paper 
tape,  loads  it,  sorts  the  data,  then  plots  the  first  pair  of 
characters  in  each  line  by  printing  disconnected  diamonds, 
and  finally  plots  corresponding  pairs  of  characters  by  means 
of  a  continuous  line  unless  there  is  no  value  present  and  the 
line  is  broken. 

The  memory  printouts  of  the  above  programs  are  given  in 
Appendix  A.  Two  paper-tape  copies  of  each  program  are  also 
kept  at  the  Radiophysics  Laboratory. 

The  rest  of  the  programs  are  written  in  BASIC  language 
and  carried  through  in  the  Dartmouth  Time  Sharing  System. 
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The  method  of  overlapping  polynomials  is  performed  essen¬ 
tially  by  two  programs:  "JET"  which  calculates  a  fixed  set  of 
coefficients  corresponding  to  the  magnetic  latitude  of  the 
Dartmouth  ionosonde  network  and  need  be  done  only  once,  and 
either  "JET-1"  or  "JET-2"  which  convert  virtual  heights  into 
true  heights  and  present  the  results  in  the  form  of  iso-height 
contours  with  different  grades  of  detail. 

The  "TID"  program  takes  a  set  of  two  contours  per  station 
and  by  means  of  crosscorrelation  and  autocorrelation  techniques, 
determines  the  main  T.I.D.  characteristics:  speed,  azimuth, 
wavefront-inclination,  power  spectra,  and  comparison  with  the 
ionospheric  predictive  function. 

The  BASIC  programs  are  also  given  in  Appendix  A  and  two 
tapes  of  each  are  kept  at  the  Radiophysics  Laboratory. 
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B„  OBSERVATIONS 

The  Dartmouth  three-station  ionosonae  network  commenced 
operation  late  in  1968.  All  of  the  1969  ionograms  have  been 
visually  inspected  by  motion-picture  projection.  Of  some  38 
likely  events#  only  7  turned  out  to  be  travelling  ionospheric 
irregularities  and  these  are  treated  in  the  next  subsection. 

The  possibility  of  a  simultaneous  observation  with  the 
Canadian  satellite-borne  ionosondes  was  actively  pursued  on  a 
suggestion  of  Prof.  Morgan.  The  actual  results  were  not  satis¬ 
factory  but  the  basic  ideas  of  much  of  the  work  described# evolved 
from  that  endeavor.  The  lack  of  success  was  essentially  due 
to  the  poor  quality  of  the  satellite  data  used:  by  the  Canadian 
laboratory's  own  standards#  most  of  the  data  qualified  as 
guesses,  their  lowest  contours  were  quite  distant  in  height  from 
our  highest  contours  so  as  to  render  the  application  of  the 
ionospheric  predictive  function  quite  doubtful,  and  finally, 
the  data,  theirs  and  ours,  were  presented  in  the  form  of  iso¬ 
ionic  contours  which  required  inversion  of  the  ionospheric  pre¬ 
dictive  function  for  comparison  with  Hooke's  theory. 

Dr.  J.H.W.  Unger  of  the  Bell  Telephone  Laboratories  called 
to  our  attention  the  fact  that  Dr.  J.V.  Evans  had  modified  the 
incoherent  backscatter  radar  at  the  Millstone  Hill  Observatory 
in  northeastern  Massachusetts  to  obtain  electron-density  height 
profiles  up  to  600  km  or  more  in  two  minutes  and  by  this  means 
had  observed  one  substantial  T.I.D.  This  provided  an  excellent 
opportunity  to  test  our  measurements  and  related  assumptions, 
and  a  comparative  analysis  was  undertaken  immediately  and  is 
reported  in  subsection  2. 
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1.  Routine 

As  has  already  been  said,  the  analysis  of  all  the  iono- 
grams  corresponding  to  the  year  1969  resulted  in  finding  seven 
events.  It  should  be  noted  that  occasionally  a  likely  T.I.D. 
had  to  be  discarded  when  one  of  the  three  stations  was  not  in 
proper  operating  condition. 

The  general  treatment  of  the  data  has  also  already  been 
described  but  a  brief  description  of  the  velocity  determination 
and  spectral  analysis  are  in  order  here. 

The  T.I.D.  is  assumed  to  be  essentially  a  wave  packet 
centered  about  a  dominant  period  component,  the  phase  velocity 
of  which  is  constant  in  magnitude  and  direction  throughout  the 
space  and  time  of  observation. 

Any  change  in  the  shape  of  the  T.I.D.  is  considered  slow 
enough  that  its  effect  on  time-delay  determination  is  negligible. 

The  techniques  used  to  determine  the  various  speeds  are 
basically  those  of  JONES  (1969)  and  are  pictoriallv  presented 
in  Figure  10.  Corresponding  pairs  of  iso-height  contours  of 
the  three  stations  arc  first  stripped  of  any  linear  trends  that 
may  be  present.  The  apparent  speeds  vA1  and  vA2,  at  a  given 
height,  are  determined  by  dividing  the  distances  from  Highgate 
Springs  and  Errol  to  Hanover  by  their  corresponding  time  lags. 

The  apparent  horizontal  speed  vAK  is  then  calculated  from  the 
geometric  construction  shown  in  the  upper  part  of  figure  10. 

The  apparent  vertical  speed  at  each  station  is  computed  by 
dividing  the  dj  f ference  in  height  of  the  iso-height  contours 
by  their  corresponding  time  lag  obtained  at  each  of  the  three 


:  35 

stations  and  an  average  of  the  three  speeds  is  taken  as  vAV, 
the  apparent  vertical  speed  of  the  T.I.D.  Finally,  the  true 
speed  v  is  obtained  from  v^H  and  by  means  of  the  geometric 
construction  depicted  in  the  lower  part  of  Figure  10. 

Enough  data  are  scaled  so  as  to  completely  encompass  the 
T.I.D.  in  time.  After  subtracting  constant  and  linear  terms, 
zeros  are  added  beyond  the  time  limits  of  the  data  in  order  to 
increase  the  recurrence  period  in  a  standard  Four:  jr  analysis 
and  obtain  finer  detail  by  closer  harmonic  spacing.  The  power 
spectra  of  two  selected  contours  at  the  Hanover  station  are 
obtained  and  the  ionospheric  predictive  function  for  the  spectral 
frequency  range  that  contains  75%  of  the  total  energy  is  com¬ 
puted.  This  experimental  predictive  function  is  then  compared 
with  its  theoretical  estimate. 

The  following  is  a  list  of  the  T.I.D.'s  recorded: 


1969  E.S.T.  (U 
5 

.T.  minus 
hrs) 

V 

4> 

T 

KP 

___ 

Jan  13,  10:00 

-  14:00 

45.6 

164.6 

74.7 

59 

Mar  02,  10:00 

-  14:00 

64.6 

34.3 

15.2 

50 

1 

Apr  29,  13:30 

-  17:30 

247.2 

171.3 

64.6 

116 

3 

Sep  05,  09:00 

-  13:00 

327.1 

169.7 

19.4 

100 

5  - 

Nov  15,  11:00 

-  15:00 

55.2 

156.6 

74.6 

190 

0  + 

Nov  17,  08:16 

-  12:16 

52.2 

179.8 

66.8 

190 

0  + 

Dec  10,  10:00 

-  14:00 

61.9 

180.5 

65.9 

67.9 

1  + 

Here  v  is  the  true  speed  in  m/s ,  <f>  is  the  East-of-North  azimuth 
in  degrees  towards  which  the  T.I.D.  was  travelling,  ij>  is  the 
inclination  of  the  wave  normal  below  the  horizontal  in  degrees, 
T  is  the  dominant  period  in  minutes  and  K  is  the  planetary 

Jr 

3-hour  magnetic  indv.  /  A  rough  estimate  of  the  accuracy  in  the 
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measurements  of  a  T.,I.D.  with  average  characteristics  yields  a  10° 
uncertainty  in  the  azimuth  and  inclination,  and  a  20%  tolerance  in 
the  speed.  The  error  involved  in  the  determination  of  the  dominant 
period  should  be  around  2.5%. 

No  statistical  analysis  of  such  a  small  population  may  be  made. 
However,  some  points  should  be  noted.  The  fact  that  most  T.I.D.'s 
occurred  at  around  noon  merely  reflects  an  observational  bias  since 
the  electron  density  variations  are  proportional  to  the  equilibrium 
density  value  and  this  achieves  a  maximum  at  about  noon  time.  The 
speeds,  inclinations,  and  dominant  periods  are  all  well  within  the 
range  of  observed  T.I.D.'s  reported  by  other  workers.  With  regard 
to  the  azimuth  of  the  irregularities  it  should  be  noted  that  all 
but  one  we re  travelling  essentially  southward,  not  in  disagreement 
with  other  reports,  see  for  example  GEORGES  (1968).  Also,  except 
for  the  T.I.D.'s  of  April  and  September,  all  events  took  place 
under  extremely  quiet  magnetic  conditions.  A  sudden  commencement 
was  recorded  on  September  5,  1969  at  00:34  E.S.T. 

Figures  11  through  59  show  the  critical  frequencies,  iso¬ 
height  contours,  power  spectra  and  ionospheric  predictive  functions. 
The  iso-height  contours  have  not  been  smoothed  and  their  finest 
roughness  corresponds  to  the  two  minute  data  sampling  period. 

As  can  be  seen  from  the  figures ,  the  critical  frequency  or 
maximum  ordinary-ray  frequency  reflected,  foF2,  varies  in  a  manner 
closely  following  the  highest  iso-height  contours.  This  is  because, 
at  times  of  large  ionization,  the  height  of  the  maximum  electron  den¬ 
sity  does  not  change  markedly  during  the  passage  of  a  T.I.D.  although 
the  magnitude  of  the  density  at  the  maximum  does  change.  Thus,  the 
foF2  curve  is  equivalent  to  an  iso-height  density  contour  at 
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the  height  of  the  maximum.  Because  the  echo  delay  sharply 
increases  as  the  ionogram  approaches  foF2,  that  frequency  is 
a  parameter  which  can  be  readily  followed  when  scanning  the 
ionograms  as  accelerated  motion  pictures,  and  periodic  vari¬ 
ation  of  it  is  a  clue  to  the  presence  of  a  T.I.D.  However, 
because  the  plasma  frequency  is  proportional  to  the  square 
root  of  the  electron  density,  a  perturbation  of  electron  density 
such  as  ten  percent,  appears  as  about  half  as  much  in  frequency. 
Furthermore,  at  night  the  unperturbed  electron  density  is  an 
order  of  magnitude  less  than  in  the  daytime  so  that  the  fre¬ 
quency  is  about  one-third  as  much  and  a  ten  percent  variation 
in  electron  density  is  but  two  percent  in  frequency. 

Before  the  advantage  of  presenting  the  data  as  iso-height 
contours  had  become  apparent,  the  seven  T.I.D.'s  found  in  the 
1969  data  had  been  reduced  and  presented  as  iso-ionic  contours. 
Because  this  is  the  conventional  method  of  presentation,  these 
results  are  included  in  Figures  50  through  101. 


2.  T. I. D.  Observed  at  Millstone  Hill 

EVANS  (1970)/  using  special  techniques  to  obtain  N(h) 
profiles  in  two  minutes  with  the  incoherent  scatter  radar  at 
the  Millstone  Hill  Observatory  in  northeastern  Massachusetts, 
has  so  far  detected  one  substantial  T.I.D.  A  particularly 
clear  section  of  this  T.I.D. ,  first  noted  by  him  and  called 
to  the  author's  attention  by  J.H.W.  Unger  of  the  Bell  Telephone 
Laboratory  (Whippany,  N.J.),  occurred  from  about  20:00  to  22:00 
E.S.T.,  May  28,  1970  and  is  the  subject  of  this  subsection. 

Due  to  the  low  values  of  electron  densities  at  the  time,  it 
is  not  feasible  to  calculate  true  heights  from  the  ionograms 
recorded  by  the  Dartmouth  ionosonde  network,  and  virtual  heights 
corresponding  to  2  and  2.5  MHz  are  used. 

Figure  102  shows  the  set  of  virtual  height  iso-ionic  con¬ 
tours  recorded  by  the  ionosonde  network  at  2  and  2.5  MHz.  The 
U.T.  interval  shown  corresponds  to  the  local  standard  time 
interval  mentioned  above.  In  Figure  103,  the  corresponding 
true  height  iso-ionic  contours  obtained  at  the  Millstone  Hill 
station  are  seen. 

The  horizontal  direction  and  apparent  horizontal  speed, 
are  obtained  from  a  cross-correlation  of  the  2  MHz  contours 
of  the  Dartmouth  stations,  while  the  apparent  vertical  speed 
is  computed  from  a  cross-correlation  of  the  two  Millstone  Hill 
contours . 

With  the  horizontal  direction  and  apparent  horizontal 
speed,  it  is  possible  to  determine  spatial  variations  of  the 


39 


T.I.D.  Figure  104 shows  the  virtual  height  iso-ionic  contours 
for  2  and  2.5  MHz  over  Hanover  at  02:00  U.T.,  May  29,  1970, 
assuming  invariance  in  the  lateral  direction  and  showing  +  500  km 
of  lateral  extent. 

The  dominant  period  of  the  T.I.D.  was  92  min;  its  true 
speed,  332  m/s;  the  azimuth  <J>  toward  which  it  was  moving,  194° 
East  of  North;  and  the  inclination  of  the  normal  to  the  wave- 
front,  29.5°  below  the  horizontal l  From  the  dominant  period 
and  apparent  horizontal  veloc  cy  which  was  381  m/s,  a  dominant 
apparent  horizontal  wavelength  of  some  2100  km  is  obtained. 

The  predicted  lag  in  the  observation  times  at  Millstone 
Hill  relative  to  Hanover  may  be  calculated  from  the  geometry 
shown  in  Figure  1Q§  and  it  is  found  to  be  5  minutes.  However, 
the  actual  lag,  found  by  cross-correlating  the  2  MHz  contours  at 
Hanover  and  Millstone  Hill,  is  10  minutes. 

After  considerable  effort  dedicated  to  explaining  this 
discrepancy,  it  seems  that  values  of  foF2  obtained  with  an 
ionosonde  at  Millstone  and  used  for  normalizing  the  radar,  are 
in  question.  Because  a  periodic  variation  appears  in  foF2  as 
a  T.I.D.  passes,  use  of  erroneous  values  can  introduce  a 
phase-shift  in  the  radar  data.  Unfortunately,  none  of  the 
rtations  of  the  Dartmouth  network  provides  unambiguous  values 
of  foF2  during  the  period  under  consideration  and  neither  does 
the  ionosonde  at  the  Lowell  Technological  Institute  near 
Millstone.  However,  MONTES  (1971),  using  a  four-station  doppler 
sounding  array  in  northern  New  Jersey,  has  found  this  T.I.D.  in 
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his  data  and  obtains  a  speed  and  azimuth  in  substantial  agree¬ 
ment  with  our  values.  Unfortunately,  due  to  the  low  electron 
density,  he  has  data  only  on  one  of  his  frequencies  and  cannot 
provide  information  with  which  the  time-lag  relative  to  Hanover 
could  be  obtained.  The  observation  of  the  event  by  his  station 
array  as  well  as  by  the  Dartmouth  stations ,  establishes  the 
lateral  extent  of  the  T.I.D.  wavefront  as  being  at  least  320  km. 

The  ratio  )N'  (z^)  j ,  determined  experimentally,  and 

its  theoretical  value  determined  from  equation  (49),  have  been 
compared  for  this  event.  Five  iso-height  contours  for  231,  240, 
249,  258,  and  267  km  recorded  by  EVANS  (1971)  with  the  Millstone 
radar,  have  been  used.  In  Figure  106,  two  iso-height  contours 
corresponding  to  240  and  258  km  are  illustrated,  while  Figure  107 
shows  their  corresponding  normalized  power  spectra. 

In  the  computation  of  equation  (49),  the  equilibrium 
density  values  were  taken  as  the  average  values  at  the  240 
and  258  Ion  contours.  The  slopes  at  240  km  and  258  km  were 
obtained  by  taking  the  difference  of  the  average  values  at 
249  and  231  km,  and  267  and  258  km,  respectively,  and  dividing 
by  18  km.  The  value  0.5  x  10-5km-1  has  been  taken  for  kz^  (=  ^~) . 
The  assumption  of  a  constant  phase  velocity,  allows  the  compu¬ 
tation  of  kj^  as  a  function  of  frequency. 

In  order  to  reduce  the  effects  of  noise,  only  a  band  of 
frequencies  where  most  (75%)  of  the  energy  of  the  signal  lies, 
is  used.  The  result  of  dividing  the  experimental  ratio 
)n'  (z^)  j  by  its  theoretical  estimate  is  shown  in  Figure  108. 
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Inspection  of  Figure  108  reveals  good  agreement  between 
theory  and  experiment  throughout  most  of  the  frequency  range. 
It  can  be  seen  that  the  ratio  of  the  theoretical  predictive 
function; as  derived  from  Hooke's  theory, to  the  experimentally 
determined  predictive  function, is  quite  close  to  unity  at  all 
points.  The  small  apparent  over-all  attenuation  is  of  little 
significance  for  it  could  be  easily  removed  by  adjusting  the 
value  of  k  . .  The  small  emphasis  at  both  ends  of  the  spectrum 

Zi> 

may  be  due  in  part  to  the  broadening  of  the  signal  spectrum, 
and  in  part,  to  residual  non-T.I.D.  variations  in  the  data. 


IV.  CONCLUSIONS 


In  the  theoretical  aspect  of  this  work,  there  have  been 

two  achievements.  Firstly,  the  characteristics  of  the  gravity 

wave  resonant  mode  ic  •  S  =  0  have  been  studied.  Its  existence 

o 

should  now  be  sought  experimentally  either  as  T.I.D.'s  or  perhaps 
as  some  permanent  feature.  In  the  case  of  non-equatorial  lati¬ 
tudes,  the  bandpass  feature  of  this  mode  could  be  detected. 
Secondly,  using  Hooke's  results,  the  concept  of  an  ionospheric 
predictive  function  is  introduced.  If  accepted  by  other  workers 
in  the  T.I.D.  field,  the  iso-height  contour  presentation  of 
T.I.D.  data  should  become  standard  instead  of  the  traditional 
iso-ionic  presentation.  Although  experimental  agreement  with 
theoretically  determined  predictive  functions  has  been  found 
to  be  very  good,  the  basic  theoretical  model  used  cou  1  be  pro¬ 
gressively  improved  by  including  the  effects  of  viscosity,  a 
constant  background  wind,  and  by  relaxing  the  adiabatic  condition. 

In  the  experimental  aspect,  the  data  processing  system  has 
been  brought  into  full  operation  and  seven  large  T.I.D.'s  have 
been  found  and  analyzed.  The  data  processing  system  could  be 
refined  further  by  changing  from  paper  to  magnetic  tape  or,  even 
better,  the  SPIRAS-65  computer  could  be  directly  linked  to  the 
Dartmouth  Time  Sharing  System,  thus  avoiding  a  few  steps.  The 
Dartmouth  Time  Sharing  System  is  used  for  the  true-height  com¬ 
putations  and  the  spectral  computations.  Now  that  the  SPIRAS 
computer  has  an  8k  memory,  the  spectral  computations  could  be 
done  with  it.  However,  it  is  unlikely  that  it  could  accommodate 
the  true-height  computations  without  still  more  memory.  A 
fourth  ionosonde  station,  possibly  semi-mobile,  would  be  highly 
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desirable.  Inasmuch  as  the  observed  T.I.D.'s  travelled  pre¬ 
dominantly  from  north  to  south,  it  could  be  placed  north  or 
south  of  the  three-station  network  to  check  the  propagation 
of  T.I.D.’s  with  theoretical  predictions  based  upon  the  param¬ 
eters  determii.ad  experimentally  by  the  three-station  network. 

It  could  also  be  placed  east  or  west  of  the  three-station 
network  to  investigate  the  lateral  extent  of  T.I.D.'s. 

The  routine  search  for  T.I.D.'s  should  be  continued  in 
order  to  obtain  a  sufficiently  large  number  of  events  for 
statistical  analysis.  Though  not  included  here,  a  BASIC  pro¬ 
gram  has  recently  been  developed  to  search  automatically  for 
a  T.I.D.,  given  the  values  of  foF2  scaled  by  means  of  the 
SPIRAS-65  program,  "critical  frequencies"  and  this  should 
greatly  assist  in  finding  nighttime  events  and  smaller  daytime 
events  than  those  discussed  here. 

The  Dartmouth  ionosonde  network  is  especially  well-suited 
for  the  continuous  monitoring  of  the  ionosphere  at  a  particular 
time  of  day  (noon,  sunrise,  sunset,  etc.).  Understanding  iono¬ 
spheric  behavior  at  these  times  may  provide  additional  insight 
into  the  T.I.D.  field. 


V.  APPENDIX:  COMPUTER  PROGRAMS 


A  total  of  eight  computer  programs  are  presented  here. 
The  first  four:  "Ionogram  Scaling",  "Plot",  "2-virtual 
heights/ionogram" ,  and  "Critical  frequencies",  are  written 
in  SPIRAS-65  machine  language  and  appear  under  the  headings: 
"SPIRAS",  "SPIRES",  "SPIRIS" ,  and  "SPIRUS" ,  respectively. 

The  leftmost  number  of  each  line  corresponds  to  the  address 
of  the  next  number  while  the  following  numbers  are  stored  in 
successive  locations.  Thus,  in  "SPIRAS",  000110  is  stored 
at  location  001000  while  000160  is  at  001007.  It  should  be 
noted  that  the  number  system  used  in  this  language  is  octal. 

The  last  four  programs  are:  "JET",  "JET-1" ,  "JET-2" ,  and 
"TID",  written  in  BASIC  language  and  require  no  additional 
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SPIRAS 


10 NOGRAM  SCALING 


001000 
001010 
001020 
001C3C 
001040 
001050 
001060 
001C70 
OCl 100 
0C1  1 10 
001 120 
001 130 
001140 
C01 150 
001160 
001 170 
001200 
001210' 
001220 
001230: 
001240 
001250: 
001260: 
001270: 
001300: 
001310: 
00I32C: 
001330: 
CO  13  40 : 
001350: 
C01360: 
C  0 1 3  7  0 : 
00 1400: 
001410- 
001420: 
001430: 
001440: 
001450: 
CO  1 460: 
001470: 
CC 1 500: 
001510: 
001520: 
001 530: 
001540: 
001  550: 
001560: 


0001 10,000222,076041 ,101 002 ,040000,177775, 1 03 1 02 » 000 1 60 
000200, 000050,000333,005140, 176176, 00CC40, 000053,004140 
176031,001041,000040,000036,026001 ,000307,000000,002161 
G01 44 5, 005 100, 176002,076025, 177773,002175,004200, 177770 
076020,076017, 177741,000171,000000,101002,000100,177775 
102102, 177771 , 002 12 7, 076005 , 076004 , 076003 , 1 77754 , 002 1 2 1 
000171 ,000000, OCl 043, C777SC, 177772,070000,072142,000000 
174000, 110000, 1 13152,050000,106152,020000,034475  044000 
136552,044000,072552,014000,  1151  12,050000,052552  6'  40  50 
01 2552, OOCOOO,OOOCOO, 000000, COCOOO,  COOOOC  ,000000,.  00(  00 
000000, 000000,000000, 000000,000000, OOOOOO.OOOOC , 1  T0000 
170512,050000,176552,000000,072142,070000,1761*  . 104000 
176552,004000, 174512, 164000,072152, 174000, 174410,000000 
107742,074000,042040, 104000,174424,100000,176040,013700 
174210,174000, 174630,O7O0CO,O72142.ClCOOO, 174512,1 70000 
072162, 1 10C0C, 174532,044000, 1 12552,004100,004176,074000 
076040,040700,03  5040, 101 70C,07603C, 052  ICC,  105210,010100 
004270,  1 1 6 1  CO,  107152,002121  ,000020, 000000, 0001 30,C00020 
101002,040100, 177775, 103102,001054,001 154, 0C7034, 167765 
027764,003774,177765, 100000,0000 10,  1  1 7757 , 1 C2 1  CO, 111154 
071044,11 1651,021777,071626,071541,001351,005200.171443 
131652,071625,071641 ,000316,001 757,003776, 177772,05 1773 
021771,1 1 1 774, C5 1772, 02 1772, 1 11776,051774,021775,1 11773 
051775,021775,002122,0002  /  l, 001776,03  1770, 111  772,002122 
000271 ,001771,031767,1 1 1 77 1 , 002 1 22 , C00270, CC22 46 , 03 1 7S6 
0001 10,000372,002 122,000271 ,001776,031765,131653, 1217S6 
013564,001401 ,023  727,003  777,  177772,071526,05 1773,021727 
002521,051774,002322,01 1770,041727,021726,002126,111726 
053  74  4,004  14  0, 1 76002 , 00600 1 , 177772, 1 13  744,053743,021353 
1 1 1726,053743,000270,000000,001352,035777, 1 13744,043743 
001 141,023744,002721 ,001042,000040,000044,001043,045777 
025777,000261 ,001777, 113744,051726,004140,176001,003777 
000240,00001 5,071626,051773,021725,002521 ,051774,021724 
002322,01 1770,041725,053727,000200,000001 ,176002, 17SC05 
177761 ,000200, 177777,177756,002000, 121725,011757,041724 
002322,01 1765,001401,004040, 1 7774 5 , 025777 ,00026 1 , OC 1 777 
003777,177740,176003,035777,000261 ,001777,0001 10,005215 
025777, C0C261 ,001777,005001 ,171521,111215,051554,005100 
171513, 1 1 121 5,001044,051555,005100, 17151 1 , 1 1 1215,001050 
051656,005100,171506,111215,001054,051233,005100,171503 
11 1215,041657, 171513,1 1 1215,041650, 171513, 11 1215,041661 
171513, 11 1215,04  1662,021215, 102  100,  1  1 1663,0  71044.071541 
1 71243 , 1  I  ICO  1 , 07 1 044 , 1 l 1 372, 071044  ,071 574, 1 1 1 44  7,0^5777 
1 1 1 777,051372,021561 , 1 1 1651 ,021777, 115777,051447,004100 
171  546,  !  1  1447,071044,001  150,07  19^1,4  ,  1  71556,  1  167  7 /.001545 
041664,071044,002121,001445,041664,071044 ,00026 1,001777 
CC01 00, COO 000, 171535, 111447, 071044,001150, 071 04 4, 0715 74 
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*?r  ;  ?<• 


1  y  <Tv» 


■  -rnypiM-'r — * 


SPIRAS  (COMT  IN’UEO  ) 


001370:  1  1 1 154,071044,000000,000171 ,000000, 131 664, 002 12 l ,071044 
001600:  003777, 177773,177770, 101005,016000, 176001 , 177774,071631 
001610:  101005,016000, 1 77770, 07 1  S3  1 , 101005,016000, 177764,101105 
001620:  103205,005100,177760,005200, 1 77756,000171,000000, 177753 
001630:  000171,000000, 1 11657,021635,000100,000000, 177775,177770 
001640:  000171,000000,000100,000000, 177775,000100,000000, 177775 
001650:  177767,002000,000006,000035,021530, 1 12600,054000,000002 
001660:  00001 0,000250, 003250, 000207, 000060,000000, 000072,0001 44 
001670:  00021 1 ,000251 ,000313,000345,000376,000431 ,000454,000501 
001 700:  000525,000551 ,000612,000654,000707,000741 ,000770,001015 
001710:  001041,001063,001104,001  123,0011-43,001  160,001176,001213 
001720:  001230,001246,000540,000740, 176140,002300, 176400,000520 
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SPIRES 


PLOT 


OOIOOO; 
001010; 
001020; 
001050; 
0C1 040: 
001050; 
CCIOSC: 
001070: 
001 100: 
CC1 110: 
001 120: 
001130: 
001 140: 
001130: 
001 i SO: 
001170: 
C012C0: 
00121 C: 
001220: 
001230: 
00124C; 
001250: 
001250: 
001270'] 
0C130C: 
OC 13 10: 
001320: 
001330: 
001340: 
001350: 
001350: 
001370: 
00140C: 
0014  10 
001420 
CO  1 4  30 
001440 
CO  1  450 
00 1  450 
OC  1  4  70 
001500 
001510 
001520 
001530 
001540 
00155 C 
C  0 1 550 


131575, 002  121, 00002 5, 001777,0 03 777, 177774*105102, 1 02 002 
CO 022 1  , COO  110, 002000,02 1777, 02 1776,002 121, 000130,000200 
022 177, 003777, 177775,071 100,051 003,005 100, 171 123,041003 
000050, 175572,005100, 171040,041031 ,022200,006001 , 171023 
1 10200,025777,000261 ,001777,000130, 177770, 12221 1,035777 
000261,001777,006001 ,177772, 171073,002000,002000,002000 
002000,002000, 002000, 002000, 002000,002000, 002000, 002000 
002000,002000,002000, 1 1 1 77S , 000040 , C00C1 4 , 1 7 1 0 13 ,000 1 7 1 
001024, 101002,040000, 177775,103102,001051,001 144,000050 
003060,101002,040000, 177775, 103202,001251,001351 ,000030 
0001 17, 041120, 171077, 102002,000023,1 1 1214,071230, 131222 
1  I  1220,07123  7, 1 1 1 2 1 5 , 071 230, 1 3 1 223 , 1 1 1 2 1 3 , 07 1237, 1 1  1 2 1 6 
021  153,  131215,1 1 1210,071237, 131221 , 1 1  121 1 ,071237,  13  1216 
1 1 1210,071237,000100,000000, 1 76001 , 1 76004, 13 1224 , 11 1212 
C71 23 7.  177757.  13  1 225 . 1 1  12 1 0. 071 23 7. 1 1 12 12 . 02 1 20 1 ,  13  12  1 5 


1  I  1212,071237, 131221 , 1 11213, 
000100,000000, 1 76001 , 176040, 
000001 ,000002,000004,000010, 
00C01 1,000012,000106,001700, 
001241,102105,101005,100000, 
071230,003777, 177775, 177772, 
071230,1 11776,051012,001550, 
002521,051213,021771, 1 11012, 
176024,051775,005100, 175014, 
004  0*0, I  76003, 1 1 1 2 i 0, C7 1 237, 

1 1 1215,071230,071335,1 1 1214, 

1 11 774,04121 3,021774, 000201, 

1  I  1266,041255,021266, 171265, 
071237, 131212,1 11216,071237, 

1  11221,071237, 131212, 1 1 1220, 
171334,131774,1 1 1212,071237, 
021373,021460,0001 1 1,004642, 
041775,021775,00272! ,005400, 
071237, 1 7^003 , 0023 54 ,111211, 
04  1255, 02 13 73, 021424,0001  1  1  , 
002324,041775,021775,071510, 
1  11213,071237,111774,041213, 
111214,071230,131774,002121, 
002002,04 12  I  0,000050,  002014 , 
1713  70, 131 775,002121  ,021  775, 

111212,071 237,000 13 0,00 14 14 , 
001434,002721 ,004 100, 176004, 
005  14  0, 176003,00233  4, 1 1 121  I , 
001545,001355,021757,031765, 
0015  75,001560,00144 1 , 04  1  757, 
000002,07123  7, 1 ! 1 2 i 3 , C7 1  23 0, 
171563,002121 , 1 7 i 542, I  1 14  24, 


071237,131216, 1 1 1212,071237 
131226, 1 1 1211 ,071237, 177757 
000020,000040,000005,000006 
000360,001 130,000310,000171 
177775,177771,000171 ,001505 
002121,021775,021774,111214 
000270, C0CO1 4 , C3 1 773 , 01 1213 
021265,0001 1 1 ,004640,005100 
002324  ,041775,021775,002  721 
176003,002334,1 11211,071237 
C 71230,131213,1 1 1213,071237 
C01 77! , 176002,176001,176031 
000171,001313,131212,11 1213 
131212, ! 1 1217,071237, 13 1212 
C71237,  131212,  1  1 1212,071237 
002121,021774,000110,002001 
CC5100, 171417,051775,002324 
171415,004040,176003, 1! 1210 
071237, 1 1 1215,071230, 1 11373 
004^42,005100, 171435,051775 
176005,1  I  1214,071230,  131213 
021774,051771 ,004100, 171417 
021774,111212,071237,000110 
005100, 175003,041463,004001 
1 1  121  I  ,071237,002000, 131222 
1 1 1210,071237,000000,000171 
1  1  1213,131213,071237,171553 
1750C! , 111210,021550,002721 
1  I  1213,021770, 1 1 1765,000114 
0C510C, 175004,002324,0001 10 
1 1 1770,0512 10, 0217 70,005100 
051255,021567,00011 1,004626 


4S 


SPIRIS 


2 -VIRTUAL  HE  1 3 HT S/I  0 NOGRAM 

OOIOOO:  0001  10,000222, 076041  , 1  01  002 , 04  0000,  1 77775 , 1 03  1 02 ,000 1  SO 
001010:  000200,000050,000333 , 005 140,1 76175,000040,000053 ,004140 
001020:  1  7*503  1 ,00104  1 , 000040 , 00003  S,  02600 1  ,  C00307 , 000035 , 002  1 S  l 
CC  l  030:  001  44  5,005100,  176002,076025, 177773,0021*75,004200, 17  7770 
001040:  076020,076017, 17774! ,000171,001507,101002,000100,177775 
001050:  102102, 177771 ,002127,076005,076004,076003, 177754,002121 
001060:  000171,001042,001043,077760, 177772,070000,072142,000000 
001070:  174000, 110000,1 13152,050000, 106152,020000,034476,044000 
001 100:  136552,044000,072552,014000, 115! 12,050000,052552,074000 
001110:  012552,000000,000000,000000,000000,000000,000000,000000 
001 120:  000000,000000,000000,000000,000000,000000,000000,170000 
001130:  170512,050000, 176552,000000,072142,070000, 176142, 104000 
001140:  176552,004000,174512,164000,072152, 174C00, 174410,000000 
001150:  107742,074000,042040, 104000,174424, 100000,176040,013700 
001  ISO:  174210,174  000, 174630,070000,072142,010000, 174512,170000 
0C1 170:  072162, 1 1 0000, 1 74 532, 04 4000, 1 12552,004100,004176,074000 
OC 1200:  076040,040700,035040, 101 700,076030,052100, 105210,010100 
001210:  C04270, 116100,107152,002121 ,000020,020000,000130,000020 
00 1220:  101002,040100, 177775,103102,001054,001154.007034,157765 
001230:  027764 , 003  774 , 177765, 100000,000010, 1 1 775  7, 1 02  I  00, 1  111  54 
001240:  071044,000110,002002,021777,071556,071545,001351,005200 
001250:  171342,0001 30, COCCCS, 07 1555,07 1545, 0003 16,001 767,003776 
001250:  177772,051773,021771 ,111774,051772,021772,11 1776,051774 
001270;  02 1 776, 1 11773,051775,021775,002122,000271,001776,031770 
CO  13 CO:  1 1 1772,002122,0002  71 ,001 77 1 , 03 1 767, 1 11771,002122,000270 
0013 10:  002246,031766,0001 10,000372,002122,000271 ,001776,031765 
001320:  131405,071556,071545,051773,021725,002521,051774,021724 
CO  1330:  121725,01 1767,041724,002322,01 1765,001401 ,000025,001777 
CC134C:  006001,171321,0001 10,000914,041777,021777,005001 ,171443 
00 1350:  105100,001550,001250, 036003,07141 1,02603 1 ,0001 10,000120 
CO  1360:  0714  H ,026024 , ^00040,000002,00 1560,000270, 000C74 ,026016 
001370:  001460,046015,001 560, 0002 70 , 00003 C, 07 | 425 , 00 1 050, 026007 
001400:  11 6005 ,071425,046004 , 1  02 1 00,  1  7 1 436 , I  77  776,000064,021000 
0C1  4  10:  0001  7!  -.001361  ,001544,001354,056005,001560,000010,000012 
00  I  420:  001 \rC , 00004 0, 000000, 177764,000171 ,001402,001560,000270 
0  Cl  *30:  COO  Cl  2, 00 124 /i  ,036001 , 000040,000 120,  177766,000!  10,000207 
CO  I  4 40:  071044,071  545, 171244, Ml COl ,07104  4,071522, 131405, l 1 1242 
CC 1 45 0:  026003, 1 1 1362,026012,000! 10,004  64  2,026004,04  1343.026014 
00  14  60:  02 ''775, 000.1  t  6, 004S2S ,  07 1  533 , 0001  CC,  CC0C1 0,  176003,07151  1 
001470:  I  1 1362.02  77  75, 1 1 i 777 , 000050 , 004 642, 004 1 00 , 177754,071511 
001500:  006001 , 177745,07151 1 ,071522,0001 10,000224,071044,000000 
COl 5  t  0:  0001 71,001503,000! 10,0002 15,071044,0001 10,000012,071044 
00 1520:  1 7776  7.000! 71,001504  ,00212! , 000 1 30 , COO  1 00, 07 1 04 4 , 003 777 
001530:  177775, 177767,00017! ,001464,001545,041613,071044,001460 
001540:  001153,041613,071044,  ”■'766,0001  71  ,001255,000100,000000 
001  550:  1777 75, 0001  00,000000, 1 777  75 , 1 7 7 76  /, COC 1 7 1 ,001 322 , 101005 
001560:  016000, 176001, i 77774 , 07 1 503 , 101005,016000, 177770,071603 


V'  TWr+lW  w  I  W*T&***ezpen*l f*»- ■(-jrp?> 
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SPI3IS  (CONTINUED ) 

001570;  101005,015000, 177764,101 105, 103205,005100, 177760,005200 
001600:  1  77756,  1 77753 ,0001  7 1 , 001 564 ,0001 10,000002,026001,000100 
001610:  000000,177775,177767,000050, 101005,016000,177764,101105 


SPIROS 


CRITICAL 

OOIOOC: 
OOIOIO: 
001020: 
001030: 
00104C: 
001050: 
0010S0: 
001070: 
001 100: 
001 1 10: 
0C1 120: 
001130: 
001  1  AC: 
001150: 
001 1  SO: 
CO  1 170: 
001200: 
CC1210: 
001220: 
CO  1230: 
C012A0: 
C0125C: 
001260: 

C 01270: 
001300: 
0C13  1C: 
CC 1 320: 
CC1330: 
001 340: 
001350: 
001360: 
001370: 
001 400: 
0014 10: 
0C1  a?0: 
001  a  JO: 
001  a<5  C: 
001  a  60* 

001400: 
001470: 
001500: 
001510: 
CC1 520: 
OC 1 630: 
0  C.  1  540: 
001550: 
001 560: 
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FREQUENCIES 


COO  1 10, 000222, 07S04 1 , 101002,040000, 1 77775 , 1 03 1 02, 0001 60 
000200,000050,000333,005140, 176176,000040,000053,004140 
176031,001041,000040,000036,026001,000307,000000,002161 
001445,005100, 176002,076025, 177773,002175,004200,177770 
076020,076017.177741,000171,000000, 101002,000100,177775 
102102, 177771 ,002127,076005,076004 ,075003, 177754,002121 
000171,000000,001 043,077750, 177772,070000,072142,000000 
174000,  1 1  COCO, 113152,050000, 106152,020000,03  44  76,044000 
136552,044000,072552,014000, 1 151 12,050000,052552,074000 

01 2552, 000000,000000, 000000,000000, 000000,000000, 000000 

OOOOCO, OCOOOC, 000000, OOCOOC, OOOCCO, OOOCCO,CCOOCO, 1 70000 
1705 12,050000, 1 76552 , COCCOC, C72 1 42 ,070C00, 176142,104000 
176552,004000,174512, 1 5400C , C72 1 52 , 1 74C00, 1 744 1 0, OCCOOO 
107749,074000,042040, 1C4000, 174424 , 1C0CC0, 176040,013700 
174210, 174000, 174430,070000,072142,010000, 174512, 1700C0 
0721 <2, 1 100CC, 174532 ,044000, 1 12552, CC4 IOC, 0C4 1 76,074000 
07604  0,04  0700,035  04  0.  101700,076030,052100.105210.010100 
004270, 1 16100, 1 07 1 52, 002 12 1, 000020, C 00000,0001 30, CCCC20 
1 01 002,04 01 00, 177775, 105102,001054,001 154,007034, 167765 
027764,003774, 177765,100000,000010, 117757,102100,111154 
071044,111651,021777,071626,07164 1 , CC 1 3 5 1 , CC52C0, 1 7 1 4 43 
13 1652, 071 <2^,07 1641 ,000316,001767,003776, 177772,051773 
021771 , 1 1 1774,051772,021772,  1  11776.051 7 74  , 02 1776, 1 1 1773 
051775,021776,002122,000271,001776,031770, 111772,002122 
000271,001771 ,031767,  11 1771  ,002122,000270,00224  6,031766 
C001 10,000379,002122 ,000271 ,001776,031765,131653,121765 
013664,001401,023727,003777, 177772,071626,051773,021727 
002521 ,051774,002322,01 1 770, 04 1 727 , 02 1 726 , 002 1 26 , 1 1 1 725 
053730,004140, 176002,006001 , 177772, 1 13730, Cc37?7 .02 1353 
1 1 1726,053727,000270,000000,001352,035777,002721 ,000050 
000015,004140, 175005,002721 , 00 1 04 4 , CCO04C , 000 i 60, 176005 
002721 ,00005  0, COCO  1 6 ,00 1 04 5 , 0COO4Q , CCC500, 04*777 , C?5 777 
000261,001777, 1 7 1 446 , COOOCO, CCOCCC, OOCOCO ,000000, COOCCO 
OOCOOC, 000000,000000, OCCOOO, COCCOC, OCOOOC, OCCOCO, OCOOOC 
OCCCO'\  000000,000000, OOOCCO, OCOCOO, COOOCO ,000000, 000000 
000000,030000, 000000,000000, OOOCCO, COCCOC, COCCOO, 000000 
CCOCCC, CCOOCC, 000000.035777,00026 1,001777,00011 0.C05215 
02^777,00026  1 , CO  1777,00500 1 ,  171521 , 1 1 ! 2  I  5 ,05 1 <54 , 005 1 00 
1  71  513,  111215 ,001  C4/i,C5  165  5,005100,  17161  1  ,  !  1  1216 ,001050 
051 656, 006  ICC, 171506, 111216,00105  4,051233,005100, 171503 
1  1  1215,04  1^57, 171513,  1  1 1215,04  1660, 171513, 1 1 1216,04 I6S1 
1 71513, 11 1215.04 1560.021215,  102 130, 111  6^3,07 1044,071 641 
l 71243, 1 i 1001 ,071 044, 1 1 1401 ,07104^,071574, 1 1 1 447,026777 
t  1  1777,06  i  401 ,02156!  ,  m*51 ,02  1777,  1  15  777,05  14  47.004100 
1  7154^,  i  i  U47,07.  ''44,001  150,071  044,  171656,  1  15777,001545 
04 16<4, 071044,002121 ,001445,04 i 5  64  ,C  7 1 C4  4 , C0C26 1 ,001777 
COO  ICO, OOCOOC, 1 71  53  5, 1114  4  7,071 04  4,001 1 50,07104  4,0715  74 


52 


SPIROS  (CONTINUED) 

001570:  ill  154,071044,000000,000171  ,000000,  131664,002121,071044 
CO  1600:  003777,177775,177770,101005,016000,176001,177774,071631 
001610:  101005,016000,177770,071631, 101005,016000,177764, 101105 
001620:  103205,005100,  177750.005200,  177756,000171,000000,177753 
001 630:  000171 ,000000, 1 11657,021635,000100,000000,177775,177770 
001640:  000171 ,000000,000100,000000, 177775,000100,000000, 177775 
001650:  177767,002000,000006,000035,021530, 112600,054000,000002 
001 SSC:  000010,000250,003250,000207,000060,000000,000072,000144 
001 670:  00021 1,000251 .000313,000345,000376,000431,000454,000501 
001700:  0005^5 , OOC55 1,000612,000654,000707,00074 1,000770,001015 
001710:  C01 041 ,001063,001 104,001 123,001 143,001 160,001 175,001213 
001720:  00 1230, 001 246, 000000,000000. 000000,000000, 000000,000000 


53 


JET 

100  ’THIS  PROGRAM  COMPUTES  THREE  SETS  OF  QUANTITIES  NECESSARY  FOR  THE 
110  ’CONVERSION  OF  VIRTUAL  HEIGHTS  INTO  TRUE  HEIGHTS  ACCORDING  TO  THE 
120  ’OVERLAPPING  POLYNOMIAL  METHOD  3Y  J.E.  TITHERIDGE. 

130  ’THE  FIRST  SET  CONSISTS  OF  TEST  FREQUENCIES  BETWEEN  2  AND  1-  MHZt 
140  ’EVERY  0.25  MHZt  AND  IS  STORED  IN  FILE  ’’FREQ". 

150  ’THE  SECOND  SET  CONSISTS  OF  AN  ARRAY  OF  FIVE  COEFFICIENTS  FOR 
160  ’EACH  FREQUENCY  AND  IS  STORED  IN  FILE  "COEF". 

170  ’THE  LAST  SET  CORRESPONDS  TO  AVERAGE  GROUP  REFRACTIVE  INDICES  FOR 
ISO  ’ALL  TEST  FREQUENCIES  AND  ASSOCIATED  PLASMA  FREQUENCIES.  THIS  SET 
190  ’IS  STORED  IN  FILE  "GRIN". 

200  FILES  FREQTsCOEFZtGRm 

210  DIM  F(49),X(5),W(5),Z(5,5),L(5,5),C(5,1)»P(49,5),A(5,1),E(49,49) 

220  FOR  1=0  TO  49 

230  LET  F(I):2+(I-l)/4 

240  NEXT  I 

250  MAT  WRITE*  1  s  F 

260  LET  HI : 1 ,4 ?4*C0S ( (90+74 )*ATN ( 1 ) /45 ) 

270  LET  H2:l ,  484*SI N ( (90+74 )*AT  N ( 1 )/45) 

280  MAT  °EAD  X 

290  DATA  .04691009,  .23076534,  .5,  .75923466,  .95308992 
3  00  MAT  READ  W 

310  DATA  .11846344,  .23931434,  .28444444,  .23931434,  .11846344 

320  DEF  FMX(F,P):Pt2/Ft2 

33C  DEF  FNT(F)  =  H2  t2/Ft2 

340  DEF  FNL(F):H1  T2/rt2 

350  DEF  FNI <F,P )= l-FNX(F.P) 

360  DEF  FMACF,P):2*FNX(F,P)*FNI(F,P) 

370  DEF  FN8  ( F,  P )  =2*F\'I  <F,P ) -FNT (F)+SQR (FNT  (F )  T2+4+FNI  (F,  P )  t2*FNL  CF)  ) 
380  DEF  FNU(F,P):SDR(1-FNA(F,P)/FNB(F,P) ) 

390  DEF  FNC(F,P):(1-FNU(F,P)  *2)  t2/FNU(F,P) 

400  DEF  FND(F,P):!/FNX(F,P)+.5*FNT(F)/FNI(F,P)T2 

410  DEF  FNE(F,P)=2*FNI(F,P) t3*FNL(F)-FNX(F,P)*FNT(F; t2 

420  DEF  FNG(F,P):FNACF,P)*FNI  CF,  P  )  +  SQR  ( FNT  CF)  t2+4*FNI  (F,P)  t2*FNL(F)> 

430  DEF  FNM (F,P ) :FNU (F, P )+FNC (F, P )* (FND (F,P )+FNE (F, P ) /FNG (F, P ) ) 

440  FOR  1=2  TO  48 
450  LET  A1  =  F(I-1)/F(I) 

440  LET  A2=F(I-2)/F(I ) 

470  FOR  J: 1  TO  3 
480  LET  F=F(I-2+J) 

490  FOR  K=1  TO  4 
500  LET  7=0 
510  FOR  R: 1  TO  5 

520  LET  T(R)=X(R)*SQR(I-(F(I-2)/F)t2) 

530  LET  G(R):F/F(I)*SQR(1-T(R)t2) 

540  LET  P:F*SQR(1-T(R)  T2) 

550  LET  7:Z  +  FNM(F,P)*T(R)*V(R)*GCR)t(K-2> 

560  NEXT  R 

570  LET  S(J,K):K*(F/F(I))t2*SQR(1-(F(I-2)/F)t2)+Z 
580  NEXT  X 
590  NEXT  J 
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JET  (CONTINUED) 

€00  M4T  2 -ZER 

€10  LET  2<l,l)=l 

€20  LET  Z(  1 ,2)=  1 

630  FOR  Mr 2  TO  5 

€40  LET  Z(M,l)r (A2- 1 )*42 t (M-2 ) 

€50  NEXT  M 

6€0  FOR  Nr?  TO  5 

670  LET  Z(N,2)r (A1-1)*A1 t(N-2> 

6  BO  NEXT  M 

690  LET  Z  (2 , 3  )  rB (  1 ,  1) 

700  LET  Z( 2 , 4) rS(2 , 1 ) 

710  LET  Z(2,5)=B(3,1) 

720  FOR  Xr3  TO  5 
730  FOR  Yr 3  TO  5 

740  LET  Z(X,Y)rR(X-2,Y-l )-B(X-2,Y-2> 

750  NEXT  Y 
760  NEXT  X 
770  MAT  LrlNV(Z) 

7<*0  LET  0(1, l)ri 
790  MAT  ArL*C 
BOO  FOR  7-1  TO  5 
BIO  LET  P(I  ,Z  )rA(Z, 1 ) 

<*20  NEXT  Z 
B30  NEXT  I 
RAO  MAT  WRITE  *2:P 
B5C  FOR  I r 4  TO  49 
B60  FOR  Jr2  TO  1-1 

B 70  LET  E(I,J)r(FNM(F(I ),F(J))+FNM(F(I),R(J-l)))/2 
BBO  IF  J<I -3  THEN  900 

B90  LET  E(I,J)r(E(I,J)+2*FNM(F(I), (F(J-l)+F(J))/2))/3 

9  00  NEXT  J 

9  10  NEXT  I 

920  MAT  WRITE  *3:  E 

930  END 


JET-1 

100  'THIS  PROGRAM  COMPUTES  PLASMA  LOG  DENSITIES  FROM  VIRTUAL  HEIGHTS 
110  'STORED  IN  FILE  ” VH4CF"  9Y  MEANS  OF  THE  TITHERIDGE  OVERLAPPING 
120  'POLYNOMIAL  METHOD.  IT  USES  THREE  SETS  OF  QUANTITIES,  PREVIOUSLY 
130  'COMPUTED  BY  "JET",  STORED  IN  FILES  "FREO",  "COEF",  AND  "GRIN  . 
140  'THE  OUTPUT  OF  THIS  PROGRAM  IS  A  SET  OF  STRINGS  EGUAL  TO  THE 
150  'NUMBER  OF  I0N03R AMS  STORED  AT  "VHJCF".  THE  FIRST  TWO  CHARACTERS 
160  'OF  EACH  STRING  REPRESENT  THE  VALUE  OF  THE  CRITICAL  FREQUENCY 
170  'AND  THE  FOLLOWING  PAIRS  REPRESENT  PLASMA  LOG  DENSITIES  FOR 
ISO  'HEIGHTS  OF  200,220,240, .. .KM 
190  MARGIN  242 
200  LET  H=1 .4S4 


2 10  FILE  #  1  * " FREQ" 

220  FILE  #  2 :"C0EF" 

230  FILE  9  3 : " GR I N" 

240  FILE  #  4t"VHYCF" 

250  DIM  F(49),?(49,5),E(49,49),A(120),V(49),H(49),K(49>,N(49),D(20) 

2  60  MAT  READ  #1»F 
270  MAT  READ  #2:P 
2 SO  VAT  READ  *3:E 

290  IF  END?4  THEN  1000 
300  LINPUT  M :  AS 
310  CHANGE  AS  TO  A 

320  LET  C=((A(l)-4S)*32-A(2)-4S)/64 

330  IF  C=0  THEN  9S0 

340  FOR  1=1  TO  49 

350  IF  C-F( I )<=0  THEN  370 

360  NEXT  I 

370  LET  N  =  I -  1 

3  SO  MAT  VrZER 

390  FOR  J=3  TO  A(0)  STEP  2 

400  LET  VC(A(0)-J+!)/2)=(ACJ^-4S)*32+A(J+l)-4S 

410  NEXT  J 

420  MAT  H=ZER 

430  FOR  1=1  TO  13 

440  LET  H(I)=V(I) 

450  "EXT  I 

460  FOR  J= 1 4  TO  N  STEP  2 
470  IF  J+ 1 > N  THEN  520 

4S0  LET  HCJ)=CVC(J+12)/2)+V((J+14)/2))/2 
490  LET  H(J+l)=V((J+14)/2) 

500  NEXT  J 

510  GO  TO  530 

520  LET  H(N)  =  V((N+14)/2) 

530  MAT  K=Z ER 
540  LET  K < 1 ) =H( 1 ) 

550  FOR  1=1  TO  4 
560  LET  R  ( I  )  =  H( I )-K ( 1 ) 

570  NEXT  I 

5^0  LET  K(2)  =  (P(2,l)*X(I)+P(2,3)*R(l)+P(2,4)*R(2)+P(2,5)'*R(3))/(l-r( 
590  LET  K (3 ) =P(3 , 1 )*X ( 1 )+P <3 ,2 )*X (2 )+? (3 ,3 )*R (2 )+P (3 , 4 )*R (3 )+P (3 , 5 )* 


1000 


;u  i\} 
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JET-1  (CONTINUED) 

600  FOR  1:4  TO  N 

610  LET  S 1 :P ( I « 1 ) *K ( I -2 ) 

620  LET  S2=P(I,2)*K(I-1) 

650  LET  C5:0 
640  LET  C4r0 
650  LET  C5:0 
660  FOR  J: 2  TO  1-2 

670  LET  C3=C3+E(I-1,J)*(K(J)-K(J-1)) 

6S0  LET  C4:C4+E(ItJ)*(K(J)-K (J-l )) 

600  LET  C5:C5+E(I+1,J)*CK(J)-K(J-1)) 

700  NEXT  J 

710  LET  S3=P(I,3)*(H(I-1)-H(1)-C3) 

720  LFT  S4:P(I«4)*(M(I)-H(1)-C4) 

730  LET  S5=P(I«5)*(H(I+1 )  —  H ( 1 )-C5) 

740  LET  K(I  ):S  1+S2+33+54+S5 
750  IF  I  <  N  THEN  770 

760  LET  K ( N) :S 1+S2+S3+S4+P ( I «  5  )*  1 , 5*H( N) 

770  NEXT  I 

7 20  FOR  1:1  TO  N 

700  LET  N(I  ):F(I > *2*1 E6/«0.5 

300  NEXT  I 

RIO  FOR  L=200  TO  2C*I NT (X ( N) /2C)  STEP  20 
S20  FOR  1:2  TO  N 
R30  IF  K ( I  )>L  THEN  350 
R4  0  NEXT  T 

S50  LrT  DCL/RO  -0)  :NC I - 1 )+ ( N (I )-N ( I- 1 ) ) /(K (I )-K (I- 1 ) )* (L-K ( I- 1 ) ) 
SSO  LET  D(L/20-S)=SC0*(L03(D(L/20-9))-L03(lCt5.5))/L0G(10t.75) 
R70  NEXT  L 

?R0  LET  H(C)=2+2*(INT(K(N)/20)-9) 

R90  FOR  1:1  TO  I  NT (K ( N) /20) -0 

900  LET  H(2*I+l)=INT(0(I)/32)+43 

910  LET  H (2* 1+2 ) : I  NT (D (I )- 32*1  NT (D  (I  )  /32 )+43) 

920  NEXT  I 

930  LET  H(l):50*C/32+4R 

940  LET  H(2):50*C-32*I NT  (50*C/32)+4B 

950  CHANGE  H  TO  BS 

960  PRINT  B3 

970  GO  TO  290 

9 BO  PRINT  "CO” 

900  GO  TO  290 
1000  PRINT  "00" 

1010  END 


try 
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JET-2 


100 
1  10 
120 
130 
140 
150 
160 

1  70 
ISO 
190 
200 
210 
220 
230 
240 
250 
260 
270 

2  SO 
290 

3  CO 
3  10 
320 
330 
340 
350 
350 
370 
3  SO 


’THIS  PROGRAM  COMPUTES  PLASMA  LOG  DENSITIES  FROM  VIRTUAL  HEIGHTS 
’STORED  IN  FILE  "VHACF"  SY  MEANS  OF  THE  T ITHERIDGE  OVERLAPPING 
’POLYNOMIAL  METHOD.  IT  USES  THREE  SETS  OF  QUANTITIES,  PREVIOUSLY 
•COMPUTED  BY  "JET”,  STORED  IN  FILES  "FREQ",  "COEF",  AND  "GRIN", 
’THE  OUTPUT  OF  THIS  PROGRAM  IS  A  SET  OF  STRINGS  EQUAL  TO  THE 
’NUMBER  OF  IONOGRAMS  STORED  AT  "VHACF".  THE  FIRST  TWO  CHARACTERS 
’OF  ^ACH  STRING  REPRESENT  THE  VALUE  OF  THE  CRITICAL  FREQUENCY 
’AND  THE  FOLLOWING  PAIRS  REPRESENT  PLASMA  LOG  DENSITIES  FOR 
’HEIGHTS  OF  HO ,H0+5 , . . . , HO+3 5  KM. 

PRINT  "HO": 

INPUT  HO 
MARGIN  242 
LET  H=1.4R4 
FILE  #  I: "FREQ" 

FILE  t  2: "COEF" 

FILE  *  3 : ” GRI N" 

FILE  #  4 :" VHYCF" 

DIM  F(49),P<49,5),E(49,49),A(120),V(49),H(49),K(49),N(49),D(20) 

MAT  READ  #1 : F 

MAT  READ  C2:P 

MAT  READ  *3 : E 

IF  E ND # 4  THEN  1020 

LINP'JT  #4:  AS 

CHANGE  AS  TO  A 

LET  C=(CA(l)-4S)*32+A(2)-4S)/64 

IF  C  =  0  THEN  1000 

FOR  1=1  TO  49 

IF  C-F( I ) <=0  THEN  390 

NEXT  I 


3  90  LET  N=  I  -  1 
400  ^AT  V=ZER 

410  FOR  J=3  TO  A (0)  STEP  2 

420  LET  V(  ( A  ( 0 )—  J+ 1  )/2)  =  (A(J)-48)*32+A(JJ-l  )  -  4  S 

430  NEXT  J 

440  MAT  H=ZER 

450  FOR  1=1  TO  13 

460  LET  H(I)  =  V(I) 

470  NEXT  I 

4R0  FOR  J= 1 4  TO  N  STEP  2 
490  IF  J+l >N  THEN  540 

500  L~T  H(J)=(V((J+12)/2)+V((J+14)/2))/2 
510  LET  H(J-*-l  )  =  V(  (J+14)/2) 

520  NEXT  J 

530  GO  TO  550 

540  LET  H (  N )  =  V (  (  N+ 1 4  )  /2 ) 

550  MAT  K=ZER 
560  LET  X ( 1 )  =  H ( l ) 

570  FOR  1=1  TO  4 
5S0  LET  R( I ) =H ( I )-K ( 1 ) 

590  NEXT  I 
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JET-2  (CONTINUED) 

500  LET  K(2>=(P(2,l>*K(l)+?(2,3)*R(l)+P<2,4)*R(2)+P<2t5)*R(3)>/<l-P<2,2)* 
610  LET  K(3)=P(3, 1 >*K( 1 )+P (3 ,2)*K <2)  +  P (3 ,3 )*R (2)  +  P (3 , 4 )*R (3 )+P (3 1 5 )*R ( 4 ) 
620  FOR  1:4  TO  N 
630  LET  S 1  :P( I » 1 )*K (I -2 ) 

540  LET  S2:P(I,2)*KCI-1) 

650  LET  C3:0 
660  LET  C4:0 
670  LET  C5:0 
6 80  FOR  J:2  TO  1-2 

690  LET  C3:C3+E(I-1  ,J  )*(X  (J)-K(J- 1)) 

700  LET  C4:C4+E(I,J)*(K(J)-X(J-m 
710  LET  C5:C  5+£  (1  + 1  ,J)*(K(J)-K(J-1)) 

720  NEXT  J 

730  LET  S3:P(It3)*(H(I-l>-H(l )-C3) 

740  LFT  S4rP(I,4)*(K(I)-H(l)-C4) 

750  LEI  S5:P(I,5)*(H(I+1)-H(l)-C5) 

760  LET  KCI  ):S1+S2+S3+S4+S5 
770  IF  I  <M  THEN'  790 

780  LET  K(N):S  l+S2+S3-*-54+P(I,5)*l  ,5*H(N) 

790  NEXT  I 

900  FOR  1:1  TO  K 

*10  LET  N(I ):F(I) *2*1 E6/«0.5 

820  NEXT  I 

830  FOR  L:HO  TO  HO+35  STEP  5 
840  FOR  1:2  TO  N 
850  IF  K ( I ) >L  THEN  870 
860  NEXT  I 

870  LET  D((L-H0)/5+l):f)(I-l>+(N<I)-K(I-l>)/(K(I)-K<I-l>)*<L-K(I-l)) 

880  LET  0 ( (L-HC) /5+ 1  ):600* (LOG  CD  C (L-HO ) /5+ l ) )-LOG( 1C  T 5.5 ) ) /LOG ( 1  Or. 75) 
890  NEXT  L 

900  LET  H(0):4+2/5*(L-H0) 

910  FOR  1:1  TO  9 

920  LET  H(2*I  H  ):INT(D(I  )/32>+43 

930  LET  H(2*I+2):INT(D(I)-32*I\’T(D(I)/32)+48) 

940  NEXT  I 

950  LET  H(l):50*C/32+48 

960  LET  H(2):50*C-32*INT(50*C/32)+43 

970  CHANGE  H  TO  3$ 

980  PRINT  3$ 

990  GO  TO  310 
1000  PRINT  "00” 

1010  GO  TO  310 
1020  PRINT  "00" 

1030  END 


100  ‘THIS  PROGRAM  COMPUTES  THE  FORMALIZED  POWER  SPECTRA  0?  A  T.I 
110  ’AT  THE  TOP  AMD  BOTTOM  OF  AN'  IONOSPHERIC  LAYER.  THIS  LAYER  I 
120  ’10  KM  THICK  AMD  IS  CENTERED  ABOUT  HO.  THE  RATIO  OF  THE  EXPE- 
130  ’R I  MENTAL  TO  THEORETICAL  IONOSPHERIC  PREDICTIVE  FUNCTIONS  IS 
140  ’ALSO  CALCULATED  USING  HOOKE’S  FORMULA.  IT  THEN  DETERMINES  THE 
150  ’SPEED  V  OF  THE  T.I.D.,  ITS  AZIMUTH  PHI,  and  us  INCLINATION 
ISO  ’PSI.  TWO  ALTERNATIVE  COMPONENT  DESCRIPTIONS  ARE  PRESENTED: 

170  ’ X , Y , Z  (EAST, NORTH, UP) ,  AND  1,2,3  WHICH  CORRESPOND  TO  EAST, 

1*0  ’MAGNETIC  FIELD  DIRECTION,  AND  THE  THIRD  DIRECTION  MAKES  THE 
IDO  ’SYSTEM  CARTESIAN.  T  IS  THE  DOMINANT  PERIOD  OF  THE  T.I.D.  AND 
2  CO  ''JA.UB.’JC  are  THE  APPARENT  VERTICAL  SPEEDS  AT  HANOVER,  HIGH- 
210  ’GATE  SPRINGS  AND  ERROL  RESPECTIVELY. 

220  VAR3IN  242 
230  FILE# l :" OPTIMO" 

24C  DIV  A(n),B(2,3S3),H(<,12l),G(2,S),"(S,l21),D(5,30) 

2  50  LET  ? 4 : A T N ( 1 ) 

9<c  PRIMT  "STARTING  COLUMN” : 

?7C  INPUT  ID 

2^0  PRINT  "AVERAGE  HEIGHT”: 

29 C  INPUT  HO 

3- DC  FOR  J=1  TO  3 <3 

3  10  LINR'JT  •  1 : At 
3  20  CHAV'S*  AS  TO  A 

33C  LET  B(1,J)=(A(2*IC+1)-4R)*32+a(2*I 0+2 )-  4G 

2  AO  LET  3(2  ,J)=(A(2*ICJ-5)-AR)*324-4(2*IC+S)-A^ 

3^0  LET  "1=?( 1 , J)/3S3 

3*C  LET  v2  =  (  (A(2*I0^3)-at)  k394-A(2*:iQ4.^)-4'7)/3S3 
37C  LET  v3=?  ( 2,J  )  /3 ?2 

3  70  LET  v4:(  (A(2*I0j-7)-4^)-*X2+A(2-*=I0j-'7)-49)/3S3 


3  -0 

:'Eyr 

j 

A  CD 

LET 

"1= 10T(5.5+.75*V1 /SCO) 

A  1C 

LET 

”2=  1 0  T  <  ^  .54..75*'!9/<oC) 

420 

LET 

"3=  1C  •  (5,5-*..75*'-3/S00) 

A30 

LET 

"i:  1 0  *  ( 5 , 5+  .  75*  \‘A  /*0C  ) 

it') 

poo 

J=1  TO  121 

A  50 

LET 

H( 1  ,J  )=°(  1  ,J  ) 

^  CC 

LET 

H  ( 9  ,  J  )  -  n  ( 9  ,  ,J  ) 

t  70 

LET 

H (3  ,  J ) =" ( l.J+121) 

flQ 

L~T 

U(A,  J)z3(2,j4-121  ) 

t*  c 

LET 

H(5,J)=9(  1  ,j4-2£9) 

500 

LET 

H(*,J)r,a(9,.j+9A9) 

5  1C 

v?yt 

J 

59C 

“OR 

1=1  TO  * 

5  3C 

Cr\o 

J=!  TO  |9| 

5  AQ 

L-  T 

S(.  ,1  )  =  G(  1 , 1  )•*•*(  I  ,.J  )/12  1 

550 

LET 

G ( 2  ,  T  ^  =  G(9,i  )+J*H  ( 1 ,  -J  )  / 1 2  ! 

5*0 

weyT 

J 

570 

»'  C  V  T 

T 

5  '-’0 

FOR 

1=1  TO  * 

590 

FO? 

y= I  TO  121 

■KflBMS-. . - 


•  CO 


ir  vvry.w  \ ;  *  -’-vp;  <. 


ll'J  (CO*.'T  I  \’UZD ) 

^00  let  7r9*(o43.3v<)/]!>0 

<i°  LET  Y=**ci EE-^EJ/lor/Ioo 

«30  "Ext  F”  »  K }  =  4<  1  .:< )"?  *  s  U  “i  )+YicQ  (? ,  I ) 

<40  «rvj  j 

^50  FOE  1=|  TO  3 
FOE  X;  -  5  T0~5 

0C0  LET  !?1“K*(!  +  SVi(K)>/2 

990  "EXT  J  --U,K  ,  )+r  (c>xl-i ,  J+'<)*F(2*I ,  J  )/J2i 

7  00  *’iYT  x 
7]0  ••'•TVT  I 
720  FOZ'lrJ  TC  3 
730  LET  Y?r-J  rj  q 
7£Q  FOE  y-Q  jj  “ 

7  50  IF  0  Cl  '< )  <vc»  TK^k!  7^o 
7^0  let  Y9  =  n(T  7  C 

7  7 C  LET  Y2r'<-r“> 

7?c  *’EyT  X 
790  LrT  X]-v9-j 
TOO  LET  X3=X2+1 
a10  LET  Y  1  =  0 (I  x J-*-5  ) 

!f2  ,L.CT  V30(I,X3-0 

L-T  ^=X1  t2*(Y9-y3)+vo  T?^fY^_v, .  ,  vr  0 

s$s  hi? 

let  uci  )r!CCO^/TCM 
770  "EXT  I  /!  <1 ) 

7C,C  FOE  1  =  4  Tq  5 

FOE  K r - J  v  TO  15 

910  LET  OCl7KllC4?-3ciK'<i f!  )J2  “ 2  |  “K*  ( 1 +S  3N' CX  > ) /2 

920  "EXT  j  CI,<  1:5>  +  r(2*I-4f J)*FC2,,J+X)/12i 

9  30  'IE, xj  x 

9^0  "ext  j 

950  FOE  1=4  TO  5 

?^0  LET  Y2r-]c1n 

9  70  FOE  X  =  o  TO  29 

9QC  IF  D(I,X)<yo  TEr"  inn 

9S0  LTT  wlw  i,'  1010 

100C  LET  V2r.x-i5 

1  0 !  0  *'Ext  x 

1020  LET  Xlrxo-I 

1030  LET  x3 :v,+! 

1040  LET  Y 1 =D (I  X I  + 1  5  } 

!c5o  i’?t  yr5^:x3+j5) 

1070  LET 

10*0  LET  TCn  =  v/5*<n  <V--Y1)  +  X.3*(YI-Y2) 

109C  MEVT  j 
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TID  (COOT  IV'JEO ) 

1  1 00  LET  L'(4)=l<3205/T(4) 

1 1 10  LET  M(5)rl47700/T(5) 

1  120  LET  ’•'=  1  .05*73 
1130  LET  SrSI'KW) 

1140  LET  CrCOS(W) 

1  1  50  LET  3:2  .****7 

1  1  EC  LET  V(<‘)  =  !\3E(,J(4)*'J(5)+S/S2R('J(4)  t2+U(5)  t2-2*  U<4  )*UC5  )*C  ) ) 
1170  LET  4:4TM(3r>R((L,(4)/U(<))+2-!)> 

1  I^C  LET  A0=(3+(1-SG*J(:J(4)))*2*?4+S3?'<U<5))*A)'M5/P4 
1  130  LET  !J(7):(’K1)*T(1  )+lK2)*T(2)+‘i(3)*T(3))/(T(l  )  +  T(2)+TC3) ) 
1200  LET  E:(,J(S)*,J(7))t2/(U(6)t2J-U(7)*2) 

1210  LET  ,J(3):Sf,R(R) 

1  22C  LET  O0:AT\’('J(£>/'J(7)  >*45/?4 
1230  LET  V(!):P/'J(E)*3I'J(4C*P4/45) 

1240  LET  V(2)=R/'J(6>*30S(A0*P4/45) 

1250  LET  V(3 )  =  -R/'J  (7) 

1 2<0  LET  I=74*P4/45 

1279  LET  V(4):VC2)*C0G(l)-V(3)*SI'i(I) 

1 2  70  LET  V(5)  =  V(2)*SIV(I)-«-V(?)*C0S(I) 

1230  OIX  V(120),YC120),P(3,121),L(242) 

1300  FOR  X:Q  TO  120 

13  10  F 0°  J: 1  TO  121-K 

1320  LET  YC-O  =  Y(K)+F(l,J)*F<ltJ+K)/90C 
1330  LET  Y(«0  =  Y(K)+r(2,J)*r  <2,J+K)/900 
134  0  ME YT  J 
1  3  50  ''EXT  K 
1 3 <0  FOR  V:Q  TO  120 
1370  FOR  X  =  0  TO  120 

13*0  LET  ?  ( t ,  Y)  :?  ( 1  .  X  ) -*-x  (X  ^CCS  ( 3*R4*  x*  X/OQC)  *2  /900 
13*0  LET  2  (*  , X ) :? (2.  X  (X  )*COG <  X*X  /*CC)*2/300 

14  00  "E  YT  •< 

14  1C  VEYT  v 

1420  "rp  v:  1  .TO  120 

1430  IF  F1>P(!  ,.v)  THE-"  14<o 

1440  LcT  =1 =?( 1 ,X) 

14  00  LET  11=1*00/* 

14^0  IF  P2>?(?,  )  THE”  14^0 
1470  LET  C2:P(2,X) 

14-7Q  LET  T2 :  1  *0C /X- 

1 4  on  v 

1  50C  FO*  Y:C  TO  120 
1510  LET  RO: PC+P ( 1 ,*) 

1520  "EYT  v 
1  530  FO2  X:0  TO  120 
1 54Q  LET  EO=ZC+?( 1 ,v) 

15  50  LET  ''0:y 

1500  IF  EO>.75-k?0  THEY  1520 
1570  -”EXT  X 

1  5 70  LET  X^:.5/(.  1  2 0+“? 4  ) 
i  03  u  LET  01  =  ((':?-"l)/10/‘' !-»-'•<*)  *2 
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TIO 


(COVTIV’JEO) 


1600 
1^10 
16?0 
1  *30 
1  *40 
1650 
1660 
1670 
1  670 
1  6<50 
1700 
1710 
1720 
1730 
1  740 
1  76C 
1  760 
1770 
17*0 
1  700 
1  *00 
1  *10 
1  *20 
1  *30 
1  *4  C 
1  *5  0 
1  *60 
1  *70 
1  **0 
1*^0 
1  9CC 
1010 
1  92  0 
1930 
1  940 
1  950 
I960 
1  97C 
1  9*0 


LET  02  -  ( ( -  V 3 )  / 1  0/\'3->-X * )  t2 
LET  S9rSIM(74*?4/45) 

OP'.  K  ( 120) 

F07  v:=0  TO  1*0 

ITT  7=V(4)  /!J(*)  ?2  *  Y  /  1  4  4  /S  9 

ITT  K  1  rS°R ( (02+ X  T2)/(D1 -*-'•<'?)) 

LET  K  (  X)  =  \'3 /'-'I* EYP (  1 0*K 9 ) *K  1 

MIT  VT  V 

FOR  vr 0  TO  120 
IF  v> v 9  THZ 1  710 

LET  ?C3,  Y)  =  X(  V)/3PR(P (2, ./)/?( 1  , X )  >*  1  CC+400 

LET  P  C 1  %:o=?(l.'*0/?l*2CO 

LET  P  ( 2 ,  X )  =  P  (2 , X  )  /?2 *2 00+200 

"TXT  Y 

LET  L  (0) =6 

J=!  TO  121 
1=1  TO  3 

L(2*I- 1 )=?(4-I, j- l ) /32+4* 
i  L(2*I)  =  P(4-I,J-l)-22*I!\’T(?(4-I,J-l)/32)  +  47 

CTVf  J 

CM 4 ^OF  L  TO  LS 
PR  I  "T  L* 

NEXT  J 
?RIVT  "CO" 

LET  V=  I  K’T  ( 1  0*  U ( * )  + .  5  )  / 1  C 

LET  A 0=  I  VT  ( 1  0*4 C4-  .6)  /  I  03^0* I  'IT  (40/36C) 

LET  ■* Or  I  "T  ( 1  C*OC+ ,  5  )  /  I  0 
FOR  1=1  TO  5 

LET  '*(  I  )  =  IVT  ( 10*V(I  )•*•  .5 )  /  I C 

"CYf  ] 

PR  I  "T  V"  Y-/3EC"  ,  "  PMI  r  ”  •  4P  .  "  Qr-r"  "  ?$  I  r  **  •  ''Q  ?  ”  DEG" 

PR  I  VT  "  VVr"  •  V(  1  )  •  ”  Y,/3EC" ,  "''Yr”  •  V(2) »"  E/SEC"  .  ”  VZ  r  "  :  V  (3  ) 1 


FOR 

FO* 

LET 
j  rr 


"  X/SE 
"  X/SE 


PRIMT  " Vlr"?V( 1 ) •" X/3EC", "  V2  =  "!V(4  )♦" X/SEC", " V3r”:  V(5) 

PR  I  V T  "  T=  "  •  I  VT  ( 5  *  ( T  1  +  T2  )  .  5 )  / 1 C  t "  XI V*  HO  r  "  •  H  0 ! "  K  X" 

FOR  1=1  TO  3 

LET  ’J  ( I )  r  I  VT  ( 'J  ( I ) *  1  c+ .  5  )  / 1  C 
ME  XT  I 

PRI  VT  "  'J4  =  "  •  *JC  1  )  ?  ”  VS  EC" ,  "  IP:"  i  U  (2 ) : "  X/S  EC" , "  UC  =  "  :  U  (3  > :  "  //SEC 
PMO 
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FIGURE  4  GRAPH  OF  l^./(l/2H)  VS  ANGLE  A 

WITH  PERIOD  T/Tg  AS  A  PARAMETER 
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FIGURE  64  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 


FIGURE  65  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 


200  ! 


CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 


CRITICAL  FREQUENCIES  AND  ISOIONiC  CONTOURS 


CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  74  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  80  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  82  CRITICAL,  frequencies  and  isoionic  contours 


FIGURE  84  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  87  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  90  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 


FIGURE  92  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 


FIGURE  93  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  100  CRITICAL  FREQUENCIES  AND  ISOIONIC  CONTOURS 
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FIGURE  10?-  VIRTUAL  HEIGHTS  AT  2  AT  2.5  MHZ 
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FIGURE  103  REAL  HEIGHTS  AT  2  AND  2.5  MHZ 

MILLSTONE  HILL 


03*00 U.T. 


CVJ 

b  o 


o  2 


o  <r> 

O  CM 


FIGURE  106  DENSITY  CONTOURS  FOR  240  AND  258  KM 
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FIGURE  1C  8  RATIO  OF  EXPERIMENTAL  TO  THEORETICAL  VALUES 
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